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Thayne Currie^, Peter Plavchan^, and Scott J. Kenyon^ 

tcurrieOcf a . harvard . edu , plavchanSipac . caltech . edu , skenyon@cf a . harvard . edu 

ABSTRACT 



We describe Spitzer IRAC and MIPS observations of the nearby 25 Myr- 

old open cluster NGC 2232. Combining these data with ROSAT All-Sky Survey 

observations, proper motions, and optical photometry/spectroscopy, we construct 

^ ! a list of highly probable cluster members. We identify 1 A-type star, HD 45435, 

with definite excess emission at 4.5-24 yum indicative of debris from terrestrial 

planet formation. We also identify 2-4 late-type stars with possible 8 yum excesses, 

Iri '. and 8 early-type stars with definite 24 /im excesses. Constraints on the dust 

luminosity and temperature suggest that the detected excesses are produced by 

debris disks. From our sample of B and A stars, stellar rotation appears correlated 

*^ ', with 24 ixm excess, a result expected if massive primordial disks evolve into 

if^ ', massive debris disks. To explore the evolution of the frequency and magnitude 

^ [ of debris around A-type stars, we combine our results with data for other young 

h>^ i clusters. The frequency of debris disks around A-type stars appears to increase 

O I from ~ 25% at 5 Myr to ~ 50-60% at 20-25 Myr. Older A-type stars have smaller 

OO 

(^ I debris disk frequencies: ~ 20% at 50-100 Myr. For these ages, the typical level 

of debris emission rises from 5-20 Myr and then declines. Because 24 /^m dust 

'k> I emission probes icy planet formation around A-type stars, our results suggest 

H I that the frequency of icy planet formation is rji > 0.5-0.6. Thus, most A-type 

stars (~ 1.5-3 M©) produce icy planets. 

Subject headings: stars: pre-main-sequence- planetary systems: formation, plan- 
etary systems: protoplanetary disks. Infrared: Stars, Galaxy: Open Clusters 
and Associations: Individual: NGC Number: NGC 2232, stars: individual (HD 
45435) 
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Introduction 



Recent studies of young stellar clusters from the Spitzer Space Telescope ( IWerner et al. 



20041 ) provide powerful insights into the emergence of mature planetary systems from pri- 
mordial disks surrounding young stars. Spitzer observations indicate that the initial stages 
of primordial disk evolution proceed rapidly. Within ~ 1-3 Myr, many optically-thick pri- 
mordial disks undergo dust settling, form inner regions cleared of dust, and undergo grain 
grow th to sizes larger than the sub-micron sizes characteristi c of the interstellar me dium 



Calvet et all boosl : llada et all bood : ISareent et"aD bood : ISouwman et aP boosh . As 



[e.g. 

these processes occur, the level of inner disk emission drops and the fre quency of primor- 



dial disks declines (e.g. [Hernandez et al.l 120071 . see also iHillenbrandl l2005r) . Primordial disks 



around early type stars disappear faster than aro und later-type stars (ICarpenter et al.ll2006 
Hernandez et al.l 120071 : ICurrie and KenyonI |2008| ) . 



By ~ 10-15 Myr, primordial disks are exceedingl y rare. The disk popu l ation is then 



dominated by optically thin, g as-poor debris disks (JHernandez et al.l l2006l : ICurrie et al. 



2007al : ICurrie and Kenyonll2008l ). Because debris dust is removed on timescales much shorter 
than the age of the star, the dust must be replenished to sustain disk emissi on. Collisions be- 
tween > km-sized planetesirnals ex cited by large-scale planet formation (e.g. iBackman and Paresce 
19931 : iKenyon and Bromleyl 120041 ) can provide this repleni shment. Though debris disks 
around early- type stars can appear as soon as ~ 3 Mv r (ICurrie and KenvonI |2008| ). the 
mid-IR debris emission does not peak until ~ 10-20 Myr (ICurrie et al.ll2008al ). This rise in 



debris emission is con sistent with planetesimal co llisions during the growth of icy protoplan- 
ets in the outer disk (JKenyon and Bromleyl l2008l ). After protopl anets reach their maximum 



size a t ~ 20-30 Myr, the level of debris emission decays as ~ t ^ (iRieke et al 



20061 ) . This decay is consistent with steady-state collision models (e.g. IWyatt et al.ll2007l ). 



12005 



Su et al. 



The epoch of maximum debris emission (~ 10-30 Myr) is also important for terrestrial 
planet formation. Planet formation models an d radiometric dating sugges t that terrestrial 
planets reach their final ma ss by ~ 10-30 Myr (IWetherill and Stewartlll993l : lYin et al.ll2002 



Kenyon and Bromleyll2006l ). Terrestrial planet formation produces debris emission observ- 
able at ~ 5-10 /im. Mid-IR Spitzer observations reveal terrestrial zone dust emission around 
many 13 Myr old stars in the massive Double Cluster, h and x Persei, and several stars in 
other 5 40 Myr old clusters/associa t ions including Sco-Cen, NGC 2547, arid the Pic Moving 



Group ( Zuckerman and Song||2004j: IChen et al 



2007bl : ICorlova et aP I2OO7I : ICurrie et al" 



2008a 



2005a 



20061: 



Lisse et al. 



200 



Currie et al.ll2007bl : iRhee et al. 



Observations of many 



^At least two other, much older stars harbor warm debri s disks that may be indicative of stochastic 
collisions in the terrestrial zone: BD +20° 307 and HD 23514 jSong et aliboost iRhee et aliboosl ) 
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warm, IRAC-excess stars in h and x Persei point to a spectral-type/stellar-mass dependent 
evolution of terrestrial zone dust, suggesting that the terrestrial planet form ation process 



runs to completion faster for high-mass stars than for intermediate-mass stars (jCurrie et al. 
2007al : ICurri]l2008h . 



Reconstructing the time history of terrestrial and icy planet formation at the 10-40 
Myr epoch requires Spitzer IRAC and MIPS observations of many clusters besides h and x 
Persei. While observations of stars in 10-20 Myr old Sco- Cen and the B Pic moving group 



provi de important studies of cold dust probed by MIPS (I Chen et al.l l2005al : iRebuU et al. 



20081 ). strong constraints on warm dust probed by < 10-15 iim broadband photometry of 
cluster members are limited to IRAS, MSX, and ground-based campaign s with much poorer 
sensitivity. The FEPS [Formation and Evolution of Planetary Systems^ iMeyer et al.ll2006l ) 
Legacy Program observed ~ 35 10-20 Myr-old Sco-Cen sources in IRAC with a range of 
spectral types (~ F8/G0 to K3). However, this sample size is too small to pro vide statistically 



robust constraints given the low overall frequ ency of warm dust emission ( iMamajek et al. 
2004J : ICurrie et al. l l2007al : ICorlova et allboOTh . 



Though h and x Persei is massive enough to investigate the frequency of warm dust from 
terrestrial planet formation in one environment, observations of other clusters are needed to 
const rain the time evolutio n of warm dust in many environments. The Spitzer survey of NGC 
2547 (JGorlova et al.l 120071 ) potentially fill this void in sampling, yielding a large population 
of cluster stars (~ 400-500) and revealing warm dust around at least ~ 6 FGKM stars. 
However, recent work has revised the age of NGC 2 547 upward from 25 Myr to ~ 35-40 Myr 
( jNaylor and Jeffries! l2006l : iMayne and Nayloiil2008l ). Spitzer studies of other ~ 15-35 Myr- 
old clusters would then clearly aid in constraining the evolution of the observable signatures 
of terrestrial planet formation. 

Understanding the evolution of cold dust emission from debris disks also requires ad- 
ditional observations of 15-35 Myr-old clusters. While 24 /xm debris disk emission around 
early-type stars peaks at ~ 10-20 Myr, the evolution of debris emission at ^ 20-30 Myr is 
not well constrained. Few clusters in this age range have been observed with Spitzer. The 
most populous cluster previously assigned to this age range, NGC 2547, is likely older. Ob- 
serving 20-30 Myr old clusters may also reveal a delay in the peak debris emission for later 
type, lower-mass stars relat ive to early-type stars as implied by planet formation models (e.g. 
Kenyon and Bromleyl 120081). Observations of y oung stars at 70 ^m may provide tentative 
evidence for this delay (JHillenbrand et al.ll2008l ). Furthermore, the disk environment around 
late type, low-mass stars may be qu alitatively different d ue to the importance of corpuscular 
wind drag in removing debris dust (jPlavchan et al.ll2005l ). Recent Spitzer studies argue that 
detectable debris disks around M stars may be rarer than around intermediate/high-mass 
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stars (JGautier et al.ll2007l : iPlavchan et al.ll2008l ). To further investigate this possibihty, MIPS 
observations of M dwarf stars provide new constraints on planet formation processes. 

In this paper, we investigate debris disk evolution and planet formation in the young 
open cluster NGC 2232 located in Gould's Belt close to the Orion Nebula Cluste r. Until very 



recently, the most complete studies of this cluster date from over 30 years ago (jClarialll972 



Levato and Malarodalll974l ). which identify NGC 2232 as a group of proper motion stars with 
B-F spe ctral types and a n age of ~ 20-25 Myr. Based on pre-main sequence evolutionary 
models. 



iLvra et al.l ( 20061) est i mate an age of 25-30 Myr and find a 25 Myr nuclear age. Both 
Clarial (119721 ) and lLyra et al.l (120061 ) show that the cluster is nearby (~ 320-360 pc). At this 



distance, IRAC observations can probe well into the M dwarf regime; MIPS observations can 
detect the photospheres of early and i ntermediate-mass stars. The cluster has an extremely 
low mean reddening (E(B-V) ~ 0.07; llvra et al.l[2006h . Thus, NGC 2232's distance, age, 
and reddening combine to make investigating debris disk evolution feasible in a way that 
complements previous studies of h and x Persei and NGC 2547. 

Our investigation is organized as follows. After describing the observations, image pro- 
cessing, and photometry of Spitzer IRAC/MIPS sources, we match sources with 2MASS/optical 
VI data in §2. Analysis of these sources shows that many may have excess emission from cir- 
cumstellar dust. To identify cluster members, we rely on ROSAT x-ray observations, proper 
motion data, optical/IR color-magnitude diagrams and spectroscopy in §3. In §4 and 5, we 
analyze Spitzer data of x-ray selected cluster members to identify stars with evidence for 
circumstellar dust, identify stars with active terrestrial planet formation, and explore any 
statistical trends in the disk population with the stellar properties and in the level of debris 
emission with time. 

NGC 2232 harbors a substantial population of stars with 24 /xm excess with a range 
of spectral types from late B to late K spectral type. At least one star (HD 45435) shows 
evidence for warm 8 //m excess. Constraining the disk luminosity and dust temperature of 
this source and comparing its SED to debris disk models shows that its warm dust emission 
is plausibly produced by terrestrial planet formation. 

This survey yields robust constraints on the evolution of cold debris disks. We identify 
a correlation between stellar rotation and 24 /zm excess for high-mass stars that may reveal 
an evolutionary link between massive primordial disks and massive debris disks. When 
combined with Spitzer data of other young clusters, our survey shows that the level of debris 
emission around A stars in the cluster is consistent with a peak in the 24 /xm debris emission 
at ~ 10-20 Myr, confirming the results of recent surveys (jCurrie et al.ll2008al ). The frequency 
of 24 /xm debris emission increases from ~ 5 Myr and peaks between 10 Myr and 30 Myr. 
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Finally, this work has important imphcations for the frequency of icy planets around 
high-mass, early-type stars. The fraction of A stars (~ 1.5-3 Mq) with 24 fim excess emission 
is high (> 50%). Because this emission likely comes from disk regions beyond the ice line 
(~ 170 K), these results imply that icy planets around A stars are common. 



2. Data 
2.1. Spitzer IRAC and MIPS photometry 

2.1.1. Observations and Image Processing 



Observations of NGC 2232 were taken with the Infrared Arr ay Camera (IRAC ; iFazio et al. 



20041 ) and Multiband Imaging Photometer for Spitzer (MIPS; iRieke et al.ll2004j ) on March 
29, 2005 and April 4, 2005 as a part of the Guaranteed Time Observations program (Program 
ID 37). Both sets of data for NGC 2232 (1= 215°, b = -7.4°) cover ~ 1.4 square degrees on 
the sky with boundaries of aaooo ~ 6'*25'"50^ - 6''30'^00" and <52ooo ~ -5°30'"00" - -4°05™00". 
The IRAC observations consist of high-dynamic range exposures of 0.6s and 10.4s at S.Gfim, 
4.5/im, 5.8/im, and 8/im. Sources were typically observed at four dither points, which yielded 
an average integration time/pixel of ~ 2.4s and 41.6s for the short and long exposures. The 
integration time/pixel was less for the edges of the region (~ 20.8 s) and more for central 
regions where dithered images overlap. The MIPS observations at 24 fim and 70 fim were 
taken in scan mode with a typical integration time/pixel of 80s. Background cirrus levels 
are low (~ 28 MJy/sr) and vary little (~ 3%) across the MIPS field. 

While the MIPS Basic Calibration Data (BCD) are largely free of artifacts, the IRAC 
BCD data require additional post-BCD processing before mosaicing. Using a custom IDL 
script written by T. C, we appli ed the appropriate array- dependent correction for point 



sources on each IRAC BCD frame (JQuijada et al.ll2004l ). Inspection of the 10.4s BCD frames 



showed that many bright stars produce severe 'striping' and column pulldown effects, espe- 
cially in the [3.6] and [4.5] channels. For columns with bright stars, column pulldown shifts 
the background flux levels by up to ~ 2 MJy/sr. We removed these artifacts by applying 
a modified version of the muxstriping and column pulldown algorithms developed by Jason 
Surace and Leonidas Moustakas available on the Spitzer Science Center website. These steps 
drastically reduce striping and bias level artifacts. Flux levels for stars change by up to ~ 
5-10% in all channels due to the array- dependent correction and up to 30% for faint stars in 
channels 1 and 2 in columns affected by bright stars. 

We then processed both the IRAC and MIPS data with the MOPEX/APEX pipehne. 
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interpolating over bad pixels. Using a bicubic interpolation with outlier rejection, the indi- 
vidual processed BCD frames were combined together into a single image in each filter for 
each exposure time. The final mosaics were inspected and cleared of any remaining image 
artifacts or pattern noise that could seriously compromise the photometry. 



2.1.2. Photometry 

Because the IRAC data are not well sampled (except for the brightest stars) and the 
point-spread function core is poorly characterize d, we performed source find ing and aperture 



photometry with a combination of SExtractor (iBertin and ArnoutsI Il996l ) and IDLPHOT 
routines in lieu of PRF fitting with APEX. For source identification with SExtractor, we 
constructed a background map for each IRAC mosaic in a 24x24 pixel area smoothed by an 
8x8 pixel box. Sources were identified by their fluxes relative to the standard deviation in the 
background rms {(Ti,kgd,rms) ■ Groups of pixels with counts greater than 5abkgd,rms {'^CFhkgd^rms) 
were identified as sources for the [3.6] and [4.5] ([5.8] and [8]) filters. The resulting source lists 
were used as input for aperture photometry with the aper.pro routine. The source flux was 
computed in a 2 (~ 2.44") and 3 pixel (~ 3.66") aperture radius, using a local background 
calculated from a 4 pixel-wide annulus surrounding each source extending from 2-6 and 3-7 
pixels. We multiplied the source fluxes by the appropriate aperture corrections from the 
IRAC data handbook, version 3.Cc|. To fine-tune the aperture correction, we compared the 
photometry for the brightest unsaturated sources derived above with that using a 10 pixel 
aperture, which should require no aperture correction. 

Photometry from the 2 and 3 pixel apertures showed excellent agreement (dispersion of 
^ 0.02-0.05 magnitudes) through ~ 13.5-14th (15th) magnitude for the [5.8] and [8] ([3.6] 
and [4.5]) channels. Beyond these limits, the dispersion in magnitudes increased to ~ 0.15 
mag. The measured pixel area of > 13th magnitude stars is completely enclosed by the 
2-pixel aperture radius. Therefore, we chose photometry from the 2-pixel aperture for stars 
fainter than 13th magnitude in all filters and the 3-pixel aperture for brighter stars. Catalogs 
were constructed for both the 0.6s and 10.4s exposures and trimmed of sources lying within 
5 pixels of the mosaic edges. 

To check the photometry, we compared aperture magnitudes with results from SExtrac- 
tor. The differences between SExtractor and IDLPHOT magnitudes are typically small, < 
0.025 mags, through 13th-14th magnitude in all filters. For the long IRAC exposures, stars 
brighter than ~ 7th-9th magnitude are saturated in all bands, with negative pixel counts in 



^http://ssc.spitzer.caltech.edu/irac/dh/iracdatahandbook3.0.pdf 
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the centers of the brightest stars. Therefore, we identified stars with m[IRAC] < 9 in any 
band from the 0.6s catalog and replaced their counterparts in the 10.4s catalog. To merge 
the final IRAC catalog with the 2MASS JHKg catalog, we employed a 2" matching radius 
between 2MASS and IRAC, resuhing in 16204, 16477, 9417, and 6087 matches in the [3.6], 
[4.5], [5.8], and [8] channels. 

The photometric errors are calculated from the Poisson error in the source counts, the 
read noise, the Poisson error in the background level, and the uncertainty in the background 
{o'bkgd,rms) ■ Figured] (top panel) shows that the source counts in the IRAC bands peak at 
[3.6] ~ 16, [4.5] ~ 16, [5.8] ~ 15, and [8] ~ 14.25-14.5. The IRAC photometry approaches 
the 10a limit at [3.6] ~ 16.5 (70.56 /iJy), [4.5] ~ 16.5 (45.13 /iJy), [5.8] ~ 14.25 (229.46 
/iJy), and [8] ~ 14 (161.09 /xJy). 

Because the MIPS point-spread function is well characterized, we performed photom- 
etry on the 24/i?Ti and 70/im MIPS data with pixel response function (PRE) fitting with 
APEX. We detect 2,702 sources at [24] and 438 sources at [70]. To identify candidate stellar 
sources, we merged the MIPS catalogs with 2MASS using matching radii comparable to 
the median position errors in MIPS: 2" for the 24 ixm channel and 4" for 70 jjnj^. There 
are 526 MIPS 24/im sources with 2MASS counterparts and 66 YOfim sources with 2MASS 
counterparts. The bottom panel of Eigure [1] shows the distribution of the 24 fiin and 70 
fim MIPS magnitudes for these sources. The source counts peak at [24] ~ 10.5 (~ 460 /iJy) 
and [70] (~ 20 mJy). The corresponding peak in the 2MASS J filter for sources detected 
at [24] is J ~ 11. The (5) 10a limits are [24] ~ (10.5) 9.75 and [70] ~ (4) 3.fl Within the 
IRAC coverage, most MIPS 24 /xm detections with 2MASS counterparts are brighter than 
J=ll but the 70fj,m data yield no 2MASS matches brighter than J=14.5. Most of the 70 
^m sources are faint in 2MASS (J ~ 16), have very red near-IR colors, and are probably 
extragalactic sources. Therefore, we restrict our analysis to the IRAC and MIPS 24 fim 
bands. Table [T] lists the full photometric catalog for sources detected in at least one IRAC 
channel or the MIPS 24 fim band. 

Many of the MIPS sources without 2MASS magnitudes are galaxies or hi ghly-reddened 



backg round stars. Using the number density of galaxies in the 24 fim filter from lPapovich et al 



( 120041 ). we expect ~ 1400-1700 extragala ctic sources in 1.4 squ are degrees. Typical fluctu- 



ations in the galaxy counts are ~ 50% (jPapovich et al.ll2004l ): thus all of the unmatched 



MIPS sources could be galaxies. Alternatively, highly-reddened background MIPS sources 



•^Setting the matching radii to the median position errors guards against false detections with the risk of 
not identifying true MIPS counterparts whose positional offsets are greater than 2" . 

"^The uncertainties quoted here do not include the zero-point uncertainty, which is « 4%. 



with 2MASS JHK^ fainter than the lOa completion hmit may lack near-IR detections. This 
region of the sky contains clusters with regions of high column density interstellar gas (e.g. 
Orion Nebula Cluster), which may have a very large distributed population as well as back- 
ground embedded clusters (e.g. Rosette Nebula). Some sources may also have positions 
determined from the MIPS mosaic that lie > 2" away from their 2MASS positions. 



2.2. ROSAT X-Ray Observations of NGC 2232 

To supplement the Spitzer data, we queried the High Energy Astrophysics Science 
Archive Research Center database (HEASAR; NASA-Goddard Space Flight Centecl) for 
ROSAT detections from the High-Resolution Imager (HRI) within 2 degrees of NGC 2232. 
Our x-ray source list is drawn from the ROSHRITOTAL, Bright and F aint Source, and Brera 



Multi-scale Wavelet catalogs (IVoges et al.lll999l : iPanzera et al.l 120031 1. Observations consist 



of two separate pointings which cover most of the Spitzer field. Exposures were taken in 
20ks, 48ks or 30ks, 36ks pairs. Two additional bright x-ray sources were observed with 0.148 
ks integrations. The HRI instrument has a resolution of ~ 5" in the center of the field. 

The ROSHRITOTAL catalog lists all sources detected by the Standard Analysis Soft- 
ware System in processed public HRI datasets. This catalog lists 299 detections, many of 
which are multiple detections of the same source. The Brera catalog consists of sources 
identified from a wavelet detection alg orithm, which accurately identifies point sources and 



extended sources ( jPanzera et al.ll2003l ). The Brera catalog contains 68 sources within the 
NGC 2232 field detected at a signal-to-noise (SNR) > 4.2. The faint and bright source 
catalogs add 35 and 5 detections. 

To find 2MASS sources with x-ray counterparts, we combined the ROSAT and 2MASS/Spitzer 
catalogs using a 5" matching radius comparable to the ROSAT positional accuracy. Merging 
the catalogs yields 79 2MASS sources with x-ray counterparts. The distribution of positional 
offsets between 2MASS and ROSAT coordinates show that most fall within ~ 2" -3" of the 
2MASS positions. In one case, a single ROSAT source falls within 5" of multiple 2MASS 
sources. We chose the source with the smaller positional offset (2.8" vs. 3.2"). The source 
3.2" from the ROSAT coordinate also has 2MASS colors indicating it is a background early- 
type star. 



^ http://heasarc.gsfc. nasa.gov/W3Browse/ 
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2.3. Optical Photometry from Iciarkl (|l972|) and iLvra et al.l (l2006h 



We also include optical UBV and BVl photometry from IClarial (Il972l ) and iLyra et al. 



( 120061 ) ■ These data 1) identify a clu ster sequence for bright stars and 2) provide a longer wave- 
length baseline. The IClarial (119721 ) photometry catalog was downloaded from the WEBDA 
open cluster database. The catalog contains ~ 43 sources with V < 12. Initial positions 
in B1950 coordinates were inspected for accuracy. Fortunately, all stars are bright enough 
to have HD numbers and precise positions from the Hip parcos and Ty cho catalogs. Based 
on these measurements, we adjusted the coordinates from IClarial (119721 ) and precessed them 
into J2000 coordinates. 



The lLyra et al.l (120061 ) survey focused on a much smaller area (~ 15'xl5' centered on the 
cluster) but wei it deep e r (to V ~ 20). The survey i dentified 1,407 stars within this region. 
We merged the IClarial (119721 ) and iLyra et al.l (120061 ) catalogs with the 2MASS and Spitzer 
catalogs, using a 2" matching radius. 



Clarial (119721 ) used multiple optical color-magnitude and color-color diagrams to identify 



stars defining a locus consistent with a cluster, to derive the inters tellar reddeni ng to each 
star, and to estimate the age of the cluster. From these methods, IClarial (119721 ) identifies 
at least 19 early-type stars consistent with cluster membership, cor nputes a me a n red dening 
of E{B — V) ^ 0.01, and estimates a nuclear age of ~ 20 Myr. iLyra et al.l (120061 ) show 
that by 'subtracting' (in V/V-I space) a field population from that containing the cluster, 
potential cluster members remain which define a clear empirical isochrone for the cluster. 
They compute an age of ~ 25-30 Myr and a slightly larger typical reddening of E(B-V) ~ 
0.07. Our motivation for including these data is to compare the colors of can didate cluster 
members in the ROSAT and Spitzer sar aples with known cluster members from I Clarial (Il972l ) 
and the locus of cluster members from iLyra et al.l (120061 ). 



2.4. Optical Spectra 

To explore the NGC 2232 stellar population in more detail, we downloaded spectra of 
stars in the NGC 2232 field from the FAST data archive housed at the Telescope Data Center 
at the Smithsonian Astrophysical Observatory. Spectra were taken by the FAST instrument 
on the 1.5m Tillinghast telescope on January 31, 2001 as a part of the FAST combo queue 
observing program (PI: Brian Patten). We also took FAST spectra of ID 18012 on April 14, 
2008 (PI: Thayne Currie). The data were taken with t he 270 g mm~^ gratin g, yielding a 
wavelength coverage of 4000-9000 A with 3 A resolution (JFabricant et al.lll998l ). We identify 
38 separate sources in the 2001 data archive, 36 of which have high enough signal-to-noise 
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to analyze. 

To derive quantitative spectral types, we first compared equivalent widths of the H I 
Bal mer lines, Call H and K, G band, Na 5890, and TiO bands with spectroscopic standards 
from lJacoby et al. I (119841 ) . In a second iteration, we computed the spe ctral types from the 
semi-automatic spectral typing code SPTCLASS (JHernandez et al.ll2004l ) and from manually 
measuring the equivalent width of indices with high correlation coefficients from smoothed 
spectra. Spectral types derived manually and with SPTCLASS agree within 0. 5-1 subclasses . 
Table[2]lists the spectral types of these stars. Many stars originally observed by lClarial (119721 ) 
were reobserved here. Agreement between spectral types derived here and those from the 
literature are good to within ~ 2-3 subclasses. 

We show two spectra in Figure [21 Both have red K^- [8] col ors but very different s p ectra l 
types. The top panel shows ID 18566, fisted as a B9.5 star in iLevato and Malarodal (119741 ) 
and classified as A2.8 here. The individual indices used to derive the spectral type (e.g. Call) 
yield estimates ranging between AO and A3. Spectral types estimated from the Balmer lines 
yield A2-A3. ID 18622, shown in the bottom panel, is a Li-poor M3.4 star with H„ emission 
(EW(IIq) ~ 8 A). Because this source has a ROSAT detection and the line is single-peaked 
and rather narrow (FWHM < 6 km s^^), the emission is likely due to chromospheric activity. 
Several other stars, all with ROSAT detections, show similar levels of Hq, emission. Thus, we 
find no strong evidence for a p opulation of late- type s tars with strong Hq, emission indicative 
of accretion (EW(H„) > lOA: Iwhite and BasrilEoOsh . 



2.5. IRAC and MIPS Colors of Sources with 2MASS counterparts 



To make a first assessment of a possible disk population among the Spitzer sources, we 
consider several color-color and color-magnitude diagrams constructed from 2MASS, IRAC, 
and MIPS data. Our goal is to isolate stars with excess IR emission from stars with pho- 
tospheric colors. For the derived age (~ 25 Myr), distance (320-360 pc), and reddening 
[Ay ^ 0.2 mag), early to mid-M type stars in NGC 2232 have J ~ 14-14.5 and intrinsic 
2MASS/IRAC/MIPS colors < 0.5. For the derived sensitivities of the IRAC data, we expect 
all cluster stars to lie above the lOcr limits for all IRAC bands. Thus, we focus our analysis 
on sources with 2MASS J < 14.5. 



NGC 2232 is a sparse cluster (jClarialll972l : lLyra et al.ll2006l ). For a normal mass function, 
the observed population of massive cluster stars impfies a total cluster population of ~ 500 
stars. Because stars with IRAC excess are rare i n young clusters with ages ~ 10-30 Myr (< 
5-10%; ICurrie et al.ll2007al : iGorlova et al.l 120071 ). we expect < 25-50 stars with IR excess. 
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There are ~ 18,000 2MASS sources with IRAC counterparts in this region. Compared 



to ma ssive clusters in relatively low background regions (e.g. h and x P^r; ICurrie et al. 



2007al ). identifying an obvious main sequence for cluster stars is challenging. Therefore, 
we concentrate on color-color and color-magnitude diagrams that allow us to isolate (i) 
background/cluster stars with a well-defined dispersion around the sequence expected for 
stellar photospheres and (ii) a handful of stars with IR excesses at least 3(j larger than the 
measured dispersion in photospheric colors. 

Figure [3] shows the J vs. Ks-[8] color-magnitude diagram and Ks-[5.8] vs. Ks-[8] color- 
color diagram (top left and top-right panels, respectively). The main distribution of Ks-[8] 
colors defines a locus that widens from 0.1 mag at J=6 to ~ 0.4 mag at J=14.5. The 
flaring of this locus to ~ 0.4 mag near the faint limit results from a combination of a) larger 
photometric errors, b) a much wider range of intrinsic colors for both cluster and background 
stars, and c) a wider range of reddening for background stars. Several sources have Ks-[8] 
colors far redder than the main distribution. There are no stars with similarly blue Kg- [8] 
colors. The main distribution of Ks-[5.8]/Ks-[8] colors extends from Ks-[5.8] = -0.1 to 0.2 
and Ks-[8] = -0.2 to 0.3 (Figure [3], top-right panel). Few sources with very red Ks-[5.8] 
colors have blue K5-[8] colors. Most sources with red Ks-[5.8] colors also have red Ks-[8] 
colors. The sources with very red Ks-[8] colors could plausibly be cluster members with 
warm circumstellar dust. 

Many sources in Figure [3] (bottom-left panel) have Ks-[24] colors greater than ~ 0.5 
and as large as 5. These excess sources lie in two groups. The first group has J ~ 6-10 with 
Ks-[24] ~ 0.5 to 2.25. If these stars are 25 Myr-old cluster members, they are massive stars 
with B0-F5 spectral types. The second group has J ~ 12-14 with Ks-[24] ~ 2-5. If these 
sources are cluster members, they are low-mass stars with K5-M6 spectral types. 

In addition to stars with IR excess, the NGC 2232 field has a clear sequence of stars with 
colors of normal main sequence stars (Figure [3l bottom left panel). This sequence begins 
at J, K,-[24] ~ 5,0 and extends to J,Ks ~ 10,-0.25 to 0.25. At J ~ 10-11, the sequence 
broadens considerably and then develops a red tail of sources with very red Kj,-[24] colors. 
With measured 5-lOcr sensitivity limits of 9.75-10.5 (§2.1), the MIPS images cannot detect 
stellar photospheres with J > 11. Thus, many K and M stars detected by MIPS are probably 
24 yum excess sources. 

Figure [3] (bottom-right panel) shows that many of the stars with photospheric Kg- [8] 
colors have excess emission at 24 ^m. The main distribution of colors ranges from Ks-[8] ~ 
-0.1 to 0.1 and Ks-[24] ~ -0.2 to 0.3. We find no evidence for sources with red Ks-[8] colors 
but blue Ks-[24] colors. About six sources have Ks-[8] > 0.3 and Kg- [24] > 1. These stars 
have substantial warm circumstellar dust, have intrinisically red photospheres (see §4.1), or 
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are heavily reddened. 

Figure E] demonstrates that some stars in the direction of NGC 2232 have robust de- 
tections of IR excess emission. The number of IR excess stars, ~ 10-20, is comparable to 
our predicted number of ~ 25-50. Thus, some of these excess stars are plausibly associated 
with the cluster. To identify bone fide cluster stars and candidate cluster stars among this 
group, we consider several methods using ROSAT detections, color-magnitude diagrams, and 
proper motions. Once we have a catalog of confirmed/candidate members, we analyze the 
disk population in §4. 



3. Identifying Cluster Members 

We now use ROSAT x-ray detections, optical/IR colors, and proper motion to identify 
candidate cluster members. The positions of ROSAT sources on optical/IR color-magnitude 
diagra ms f ollow a locus con sistent with previously identified cluster members from IClaria 
( 119721 ) and lLyra et al.l (120061 ). This locus allows us to identify additional candidate members 
from optical/IR color-magnitude diagrams and from proper motion data. 



3.1. Catalog of Cluster Members ba sed on X-Ray Active Stars and IClaria 

( |l972h 



To isolate cluster members from background s tars, we first consider the RO SAT detec- 
tions. Young late-type stars are X-ray active (e.g. iPreibisch and Feigelsonl |2005| ) . For stars 
with ages ~ 10-15 Myr, the typical X-ray flux is ^ 100-1000 times larger than fluxes for 
stars of simil ar mass with ages > 100 Myr (Chandra Observations of Orion, the P leiades, 
and h Persei; IPreibisch and Feigelsonl l2005l : IStauffer et al.ll2005l : ICurrie et al.ll2008d ). 



Figure m shows the ROSAT detections and the Claria (1972) cluster members on a J/J- 
Ks color-magnitude diagram. These stars define a clear locus consistent with the isochrone 
expected for a 25 Myr-old cluster. Bright cluster members from Claria (1972) trace the 
isochrone well for J = 6-10. For J > 10, the ROSAT detections clearly indicate the location 
of the cluster, isolating it from brighter (fainter) foreground red dwarfs (background red 
giants) with J ^ 10 and J-K^ > 0.6 and the large background population with J > 12 and 
J-K, ~ 0-0.6. 

ROSAT sources outside the 25 Myr isochrone are mostly faint background stars. Stars 
with J = 14-14.5 and J-K^ = 0.4-0.8 have bluer J-K^ colors than cluster stars at the same 
J. The apparent magnitudes of these stars are consistent with young stars at distances ~ 
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500-1000 pc and thus might be members of more distant clusters or star-forming regions. 



3.2. Candidate Members Identified from iLvra et al.l ( 120061 ;) and 2MASS 

photometry 

The ROSAT detections suggest that we can use optical or optical/IR color-magnitude 
diagrams to identify candidate members. To test th is poss i bility we compare the V/V-J 
colors for cluster members identified by ROSAT and IClarial (Il972l ) with othe r stars in the 
field. We use V-band photometry from lLyra et al.l (120061 ) and photometry from lClarial (Il972l ) 
for sources outside Lyra's coverage. 

Figure [5] shows that the ROSAT and Claria stars define a tight locus in V/V-J space. 
From these stars, we define an empirical isochrone for the cluster (solid line) with an up- 
per limit to include binaries and a lower limit to allow for slight differences in redden- 
ing/distance/age. This empirical isochrone is identical to the ISiess et al.l (120001 ) theoretical 
isochrone for 25 Myr-old > 0.5 M0 (MO-Ml) stars, brighter (bluer) than V ~ 16 (V-J ~ 
3.2). To match the empirical locus of x-ray active starqj, we shift the isochrone redward by 
~ 0.5 magnitudes for V > 16. 

Compared to Figure HJ Figure [5] has many fewer bright red stars lying above the 
isochrone, which are likely foreground late-type stars or background giant s. This differ- 
ence probably results from the smaller spatial coverage of iLyra et al.l (120061 ) . w ho focus on 



the cl uster center. To test this possibility, we compare the relative areas of the iLyra et al. 
J2006h and the IRAC /MIPS surveys. The IRACS/MIPS coverage is ~ 22.5 times that of 



the iLyra et al.l (120061 ) data. Thus, we expect 22.5 times as many bright, red stars in Figure 
Hlas in Figure [5l The observed ratio of ~ 125/4 ^ 31 is reasonably close to our expectation. 

Because cluster members define a narrow locus in V/V-J space, we identify candidate 
cluster members lying along this locus between the upper and lower bounds (-1- 0.5, -0.75 
mags; dotted lines) in Figure El We also identify stars outside the lLyra et al.l (120061 ) coverage 
but overlapping the area enclosed by confirmed members with positions on the J/J-H and 
J/J-Ks color- magnitude diagrams consistent with cluster membership (Figure [6]). Candidate 
members are drawn from a more narrowly defined locus (± 0.5 mags). We make additional 
passes through J/J-K^ and J/J-H color space and reject all candidate members that clearly 
lie outside the plausible color-magnitude range in either diagram. 



^This shift allows us to match the luminosity and colors of stars with masses < 0.5 -0.6 M(t, and to acco unt 
for uncertainties in low-temperature opacities that affect stellar atmosphere models ( Baraffe et al.lll998l l. 
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3.3. Additional Selection Criteria and the Final Membership List 

To provide a final check on our membership catalog, we compare the proper motions 
of previously identified members and other stars on the field. The median proper motion 



of the cluster is derived from Tycho observations of the early-type members in the IClaria 



( 119721 ) catalog. These stars have a mean proper motion of Ara,dec ~ -5.33 ± 1.27, -2.68 



±1.34 mas yr ^. We identify all stars within la (of their measurement err ors) as having 



a proper motion consistent with cluster membership. All members in the IClarial (119721 ) 
catalog meet this criteria. Using the Tycho catalog. Guide Star Catalog, and Henry Draper 
catalog, we identify field stars and potential cluster members. We reject three bright stars 
(IDs 2162, 15449, and 18442) and add two stars (IDs 5494 and 5625) whose proper motions 
are within la of the mean cluster motion. We also add one source from our spectroscopic 
catalog, ID 10258, which has J/J-K^ and J/J-H positions barely outside our photometric 
membership criteria. This star has a J magnitude and colors very similar to other K5-M0 
cluster members. 



In total, we identify 11 members from the IClarial (119721 ) catalog, 37 members from x- 
ray activity, and two members from proper motion. To this sample we add, 15 candidate 
members from V/V-J colors and 176 candidate members from J/J-K^ and J/J-H diagrams, 
bringing the total to 239 confirmed/candidate members. Of these candidates, 209 have 
IRAC and/or MIPS photometry. Thus, the cleaned list of x-ray active and previously iden- 
tified members are photometrically (in two color-magnitude diagrams) and, in some cases, 
astrometrically confirmed as members. The list of candidate members are all photometri- 
cally consistent with membership. Of confirmed/candidate cluster members, 32 have optical 
spectra listed in Table [21 Table [3] lists all the confirmed/candidate members with Spitzer 
photometry. 

The candidate members identified from color-magnitude diagrams contain a potentially 
large population of field stars, especially for stars with near-IR colors consistent with fore- 
ground M stars and background M giants (J-H ~ 0.6-0.7, J-K^ ~ 0.75-0.85). To better 
constrain their status, these sources need additional photometric and spectroscopic data. 
In the analysis section, we retain these sources with the qualifier that some are probably 
not cluster members. Thus, though candidate members identified from optical and near- 
IR color-magnitude diagrams may include NGC 2232 sources with IR excess, including all 
candidate members identified from these diagrams necessarily introduces uncertainties in 
determining the frequency of disks for stars later than ~ FO-GO. The samples of confirmed 
and candidate cluster members are, in order from most to least robust, the x-ray active 



sources and the IClarial (119721 ) catalog, sources selected based on V/V-J and proper motion. 



and sources selected based on J/J-K^. 
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4. Analysis of the NGC 2232 Disk Population 

4.1. IRAC/MIPS colors of Cluster Members 

We now analyze the IRAC and MIPS colors of confirmed and candidate NGC 2232 
members. Figures [7] and [8] show the Ks-[4.5], [5.8], [8], and [24] colors for these stars. 
Overplotted are the loci for photospheric colors of (pre) main sequence stars (vertical dashed 
line) and the giant locus (vertical dot dashed line); J-K^ colors for A0-M5 stars are identified 
by horizontal dotted lines. The loci for photospheric colors were derived from the STAR-PET 
interactive tool on the Spitzer Science Center website, which co mputes 2MASS-Spitzer colors 



based on the Kurucz-Lejeune stellar atmosphere models (e.g iKuruczl Il993l : iLejeune et al. 



19971 ) convolved with the IRAC and MIPS filter responses (Table jH). We also overplot th e 



reddening vector assu ming the 2MAS S/IRAC reddening laws from llndebetouw et al.l (120051 ): 
A24 is estimated from iMathis I (1l990l ) . 



For stars with photospheric near-IR emission and identical extinction, the J-K^ color sep- 
arates stars by spectral type. Cluster stars follow a well-defined locus in the near-IR/IRAC 
color-color diagrams (Figure [7]). From early to late- type stars, the full dispersion about the 
photospheric locus ranges from ~ 0.05 to 0.1 mag, roughly twice the maximum photomet- 
ric error at [5.8, 8]. There are no cluster stars with red giant colors (J-K^ > 0.9) or with 
the colors of the lowest-mass main sequence stars (> M4 spectral type). The lack of red 
giant stars suggests our selection criteria are robust. The lack of stars later than M4-M5 is 
consistent with our J-band magnitude cutoff for IRAC sources with 10 a detections. 

NGC 2232 contains a small population of stars with 5.8 fim and/or 8 fim excesses. One 
star (ID 18566) with a blue J-K^ color has a Ks-[4.5, 5.8, 8] excess (~ 0.2, 0.5, 1 mags) 
clearly larger than the zero Ks-[4.5, 5.8, 8] color expected for early-type stellar photospheres. 
The stellar spectrum confirms the A spectral type suggested by the J-K^ colors and shows 
no evidence for accretion (Figure [2]). Visual inspection of the source on the IRAC and MIPS 
frames shows no evidence for source blending (Figure [H left panel). 

While most late-type stars in the cluster have photospheric colors, a few may have K^- 
[5.8, 8] colors > 0.2-0.3 mag redder than the photosphere. IDs 6540, 9220, 18601, and 18622 
have 8 fim excesses that are > 3cr[5.8],[8] away from the right-hand vertical dotted line showing 
the 0.1 magnitude bound for photospheric sources. Thus, these stars have red Ks-[IRAC] 
colors that may indicate the presence of warm, circumstellar dust. 

One star (ID 6540) with red J-K^ colors appears to have a ~ 0.7 magnitude excess at 
8 /xm. ID 6540 may also have a borderline (~ 0.27 mag) excess at 5.8 fim. Although the 
image mosaic at 5.8 /xm does not show contamination from another source, the 8 /zm image 
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is slightly extended and may be blended. Blending could arise from the superposition of ID 
6540 with a background galaxy. Galaxies have much redder [5.8]-[8] colors than stars; thus 
a background galaxy can raise the 8 fim flux of a stellar source and remain undetected at 
shorter wavelengths. Thus, the excess for this source is questionable. 

The other late-type stars with 8 fim excess are not artificially brightened by source 
blending (e.g. ID 18601; see Figure [H right panel). Another source with very red Ks-[5.8, 8] 
colors, ID 18622, also lacks evidence for blending. The IRAC mosaics show that these stars 
have no neighbors whose light partially falls within their aperture radii. These sources are 
also ~ 12th-13th magnitude, corresponding to an IRAC flux ~ 4-15 times brighter than the 
10a limits in these filters. Thus, NGC 2232 may contain a late-type population of stars with 
warm dust emission. In contrast, all sources with J-K^ ~ 0.2-0.7 (plausibly F0-K5 stars) 
have photospheric colors. 

The population of potential 24 fim excess sources is larger (Figure [8]). Many cluster stars 
with blue J-K^ colors have 24 fim fluxes up to ~ 3-10 x their predicted photospheric fluxes 
(Ks-[24] ~ 1- 2.5). Other members with slightly redder colors (J-K^ ~ 0.2-0.5; plausibly FO- 
KO stars) may also have weak excess emission (~ 0.25-1 mags). Figure [S] (left panel) reveals 
a third population of stars with Ks-[24] ~ 0.5-1.5 that are plausibly K/M stars. While the 
errors in [24] are larger at a given magnitude than in the IRAC bands, the excesses are 
typically much larger than IRAC excesses. In §2 we showed that the number counts of 
MlPS-detected stars peak at [24] ~ 10.5 and J ~ 11. Assuming a distance of 340 pc, an 
extinction of E(B-V) =0.07, and a cluster age of 25 Myr, 11th magnitude cluster stars should 
have a G2 spectral type. Thus, the MIPS survey is likely complete to early G stars. While 
MIPS does not detect the photospheres of stars later than ~ G2, it can detect late-type stars 
with 24 fim excesses from circumstellar dust. 

All but one of the MlPS-excess sources appear to have photospheric IRAC colors (Figure 
m top-right panel). The conflrmed, early and intermediate-type (B-G) NGC 2232 members 
have a range of 24 fim excesses ranging from ~ 0.3 to 2.2 mag. However, all but one have 
a very narrow range of Ks-[8] colors (~ -0.1 to 0.1) consistent with stellar photospheres. 
The later-type sources with ~ 0.3-0.9 mag 24 fim excesses from Figure |H] (top-left panel 
and bottom-left panel) also have a narrow range of Ks-[8] colors (~ 0-0.15) consistent with 
late-type photospheres. Sources lacking 8/im excess also have a range of K5-[5.8] and Ks-[4.5] 
colors consistent with photospheres. 
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4.2. SED Modeling of HD 45435 (ID 18566): Evidence for Debris Emission 

Produced Prom Terrestrial Planet Formation 



NGC 2232 has many good candidates for debris disks. In addition to a handful of 
stars with modest IRAC excesses, there are many cluster members with 24 /xm excess 
emission. The level of this excess emi ssion (^ 0.5-2.25 mags) is much smaller than ex- 



cesses in primordia l disks (~ 5 mags; iKenyon and Hartmann 



1995: Hernandez et al. 2007 



Currie and Kenyoru 120081 ) ; however, the levels are similar to the excesses observed in de 
bris disks (~ 0.5-3 mags; iRieke et al.ll2005l : ISu et al.ll2006l : iGorlova et al.ll2007l : ICurrie et al. 



20084 



To test whether the excess sources are primordial disks, evolved primordial/transition 
disks, or debris disks, we rely on analyses of spectral energy distributions (SEDs) and the 
ratio of the observed fluxes to predicted photospheric fluxes at 8 fim and at 24 /im. Only one 
star, ID 18566 with an A3 spectral type, has clear excesses in the IRAC and MIPS bands. 
These excesses allow us to measure the temperature and luminosity of the dust. As shown 
below, our analysis demonstrates that this star harbors a debris disk. Other sources with 
IRAC (MIPS) excesses have no obvious MIPS (IRAC) excesses. Thus, we cannot analyze the 
SED in detail. Because the excesses for these sources are smaller than those for ID 18566, 
these stars probably also harbor debris disks. In §4.3, we show that the observed excesses 
and upper limits at other wavelengths are consistent with this conclusion. 



To analyze t h e SEP for HP 45435, we follow the procedure from lAugereau et al.l (119991 ) 
and I Currie et al.l (l2007bl . l2008al ). The star+disk emission is flrst flt to a sum of the stellar 
blackbody plus two scaled blackbodies from the disk to derive a) the amount of disk emission 
relative to the star and b) the temperature, and thus location, of the disk emission. We as- 
sume uncertainties of 5% (10%) in the IRAC (MIPS) bands due to possible systematic effects 
(e.g. uncertainty in the zero-point flux, phase- dependent response, and photometric errors). 
Evolved primordial dis ks such as transition disks in IC 348 and Taurus have fraction al lumi- 
nosities, Lrf/L^, KiO.l (ICalvet et al.ll2005l : ILow et al.ll2005l : ICurrie and Kenyonll2008l ). Debris 
disk s have disk luminosi t ies > 100 times smaller than evolved primordial disks ( ~ 10~^-10~^, 



e.g. IPadgett et al.ll2006l : ISu et al.ll2006l : iRhee et al.ll2007al : ICurrie and Kenvonll2008f ). Thus, 
the disk luminosity discriminates between disks of different evolutionary states. 

Second, the source SED is compa red to more rigorous d i sk mo dels. We consider 
an evolved primordial disk model from iKenvon and HartmannI (119871 ). a terrestrial zone 
debris disk model fro r a iKenvon and Bromleyl (120041 ). and a cold debris disk model from 
Kenyon and Bromleyl (120081 ). For the evolved primordial disk model, the disk tempera- 
ture scales as T oc r^^/*^. We vary the inner hole size to match the level of IRAC emis- 
sion. Debris from terrestrial planet formation is modeled at 1.5-7.5 AU from a 2.0 Mq star 
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( iKenyon and Bronileyll2004l ). For the cold de bris disk model, debris from planet formation 



is modeled at 30-150 AU from a 2.0 Mq star (JKenyon and Bromleyll2008l ) 



The min(x^) scaled blackbody fit for HD 45435 shown in Figure [TD] demonstrates that 
the amplitude of disk emission (L^/L^ ~ 5xl0~^) is comparable to the most lurninous d ebris 



disks ( e.g. HD 113766A and h and x Per-S5, see lLisse et al.ll2008l : ICurrie et al.ll2008af ) that 
the disk has very hot dust. The fits require dust populations at two different temperatures, 
~ 360 K and 815 K. This emission is consistent with warm dust in regions comparable to 
the orbits of Mercury and Earth in the Solar System. Disk emission clearly emerges above 
the stellar photosphere by 3.6 ^m and dominates over the photosphere longwards of 8 /xm. 

Comparing the source SED to predictions from transition and debris disk models con- 
firms that the HD 45435 likely harbors a terrestrial zone debris disk. The transition disk 
model drastically overpredicts the mid-IR fiux longwards of ~ 4-5 /xm. The cold debris disk 
fares better in reproducing the MIPS excess. However, it predicts only marginal 5.8 /xm and 
8 iim excesses in contrast to the strong excesses observed. The terrestrial zone debris disk 
models accurately reproduces the IRAC and MIPS fiuxes. 

In summary, HD 45435 (ID 18566) has the strongest 24 ixm excess of any NGC 2232 
member and is the only source with a strong IRAC excess. Compared to other cluster 
members, its disk is the most likely to be at an early evolutionary state, perhaps an evolved 
primordial disk/transition disk. However, SED modeling suggests that this source is instead 
a warm debris disk. The SEDs of sources with weaker 24 ^m emission are even more optically 
thin and more likely to be debris disks. Thus, optical spectroscopy and SED modeling suggest 
that the mid-IR excesses from early-type stars in NGC 2232 are due to dust in circumstellar 
debris disks. 



4.3. Frequency of 8 /i?Ti and 24 ^m Emission from Debris Disks 

We now consider the frequency and nature of the stars with excess IRAC (MIPS) emis- 
sion but no statistically significant MIPS (IRAC) excess. After defining robust criteria for 
these sources, we measure the frequency of disk emission in both bands. To place initial 
constraints on the disk properties of these sources, we use IRAC and MIPS upper limits to 
constrain the dust temperature and to derive a lower limit to the dust luminosity. These 
results demonstrate that NGC 2232 members include at least two stars with warm debris 
disks and many stars with cold debris disks. 
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4-3.1. Frequency of 8 fim Emission 

With the exception of HD 45435, NGC 2232 lacks any stars with very luminous 8 
fim excess emission (e.g. greater than twice the photospheric flux). Warm dust around 
other stars must be far weaker. The preceding section (§4.1) identifies several stars with 
Kg- [8] colors redder than 3a away from theoretical photospheric colors. However, more 
robust measures for establishing excess sources are needed. Specifically, the Kurucz-Lejeune 
synthetic Ks-[8] colors used to determine photospheric colors (and thus excess sources) are 
prone to systematic uncertainties in the model atmospheres. Young M stars, which comprise 
the majority of the 8 /im excess candidates, are chromospherically active and can have large 
fractions of their sur faces covered \yith s tarspots, causing changes in the brightness of colors 



of ~ 0.1 magnitudes (iBouvier et al.lll993l ). Thus, sources with near IR and mid-IR data taken 
at different epochs may have an intrinsically larger dispersion in colors at a given spectral 
type. This dispersion results in a larger uncertainty in measuring the amount of IR excess 
emission. Because NGC 2232 was not observed in K^ and [8] simultaneously, starspots may 
also contribute small uncertainties in identifying excess sources based on the K5-[8] color. 

To minimize the impact of systematic errors in identifying stars with 8 nm excess emis- 
sion, we determine the locus of photospheric colors empirically from IRAC fluxes. Following 



Carpenter et al.l (120061 ). we use the [4.5] filter to define the stellar photosphere and compare 
its flux to the 8/im flux to identify excess sources: log(F[8]/F[4.5]). Using the 4.5 /im flux 
has several advantages over using Kg. The scatter in the [4.5]- [8] color as a function of [8] for 
all IRAC-detected stars (not shown) is slightly smaller than that in Ks-[8]. All sources with 
8 iim detections should also have detections at 4.5 /xm. Furthermore, the theoretical Ks-[8] 
colors from the Kurucz-Lejeune stellar atmosphere models suggest that the [8] to K^ flux 
ratios as a funct ion of spectral type ca nnot be fit by a linear function, whereas the [8] to [4.5] 



flux ratios can (jCarpenter et al.ll2006l ). Finally, the 8 /im to 4.5 fim flux ratio is insensitive 



to variations in reddeni ng because stars are reddened by the same amount in either filter 



fllndebetouw et all 120051 ). 



Using 4cr iterative clipping, we calculate a best-fit line using a standard linear least- 
squares fit to the locus of cluster stars in J-K^ vs. log(F[8]/F[4.5]) space. The final linear fit 
has a slope of 0.046 and a y-intercept of -0.465 with rms residuals of ~ 3.92%. To determine 
an emprical excess criterion, we identify the threshold, in arms, beyond which there are no 
sources with 'negative excesses' (deficient flux ratios): ~ 3.4:2arms- Assuming a dispersion 
in IRAC colors dominated by photon noise, we expect less than 1 star out of 200 cluster 
members detected at 8 /xm to be more than S-iarms away from the photospheric locus. We 
adopt a 4:arms threshold to identify cluster stars with 8 /im excesses. Similar thresholds have 
been adopted by other authors to identify stars with weak excesses in both IRAC and MIPS 
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e.g. ICarpenter et al.ll2006l : iHillenbrand et al.ll2008l ) 



NGC 2232 includes four stars that are > 4crrms away from the photospheric locus in 
log(F[8]/F[4.5]) (Figure [TT]). For confirmation of these excesses, we compare their Ks-[8] 
colors to photospheric predictions. ID 18613, a cluster member identified by x-ray activity, 
has a 4.35(7 excess in the 8fim to 4.5 fim flux ratio but has a rather blue K5-[8] color (~ 
0.183) consistent with photospheric predictions. However, three stars - IDs 6540, 9220, and 
18566 - have flux ratio excesses and red Ks-[8] colors lying outside the locus of predicted 
photospheric colors. The apparent excesses for these stars are more than 3a larger than their 
photometric errors. Even though ID 6540 has a > 8cr excess, the previous section cites source 
blending as a possible source for its excess emission. In contrast, ID 18566 (HD 45435) is 
not blended and clearly has excess emission from 4.5 fim onwards. ID 9220 also shows no 
evidence for blending and has red Ks-[8] and Ks-[5.8] colors suggestive of mid-IR excess in 
addition to its excess in log(F[8]/F[4.5]). ID 6540 is identified as a candidate member based 
solely on near-IR color-magnitude diagrams; ID 9220 is identified as a member from both 
optical/near-IR and near-IR color-magnitude diagrams. Therefore, at least early-type star 
and likely one late-type star (ID 9220) in NGC 2232 have bona fide IRAC excesses indicative 
of warm dust; two other late-type stars may also have warm dust emission. Properties of 
the four 8 fim excess sources are listed in Table |5l 

While our method can identify IR-excess sources without some systematic uncertainties, 
it may conceal some bona fide excess sources with very warm dust (T^; ~ 350-600 K). Using 
the 8 fim to 4.5 /im flux ratios to define excess sources best identifies warm dust whose flux 
peaks at A > 8 fim. In terms of their photometric errors, IDs 18601 and 18622 have both 
Ks-[8] and Ks-[5.8] colors (~ 0.4 mags) greater than 3a away from the locus of predicted 
photospheric colors. However, their 8 fim to 4.5 fim flux ratios are less than 3a away from 
the empirically-defined photospheric locus. If these systems have weakly-emitting warm dust 
that peaks in the IRAC bands, then their 8 fim to 4.5 fim, flux ratios will be bluer than that 
for dust whose 8 fim, emission is on the Wien tail of a slightly colder blackbody. 

We adopt a simple model to show that the weak 8 fim, excesses around late-type stars 
yield disk luminosities consistent with debris disks. To estimate the disk luminosity at 8 
fim,, we calculate the flux of the 8 fim, excess, assume that the disk emission peaks at 8 fim, 
and originate s from a 36 K bl ackbody, and assume the stellar luminosity of a 25 Myr-old 



MO star from lSiess et al.l (J2000l ). We derive L^/L^ ra 10~'*-10~^, consistent with debris disk 
luminosities. The maximum disk luminosity from a 120 K blackbody consistent with the 24 
fim, upper limits is also L^/L^ ^ 10~^. 

To estimate the frequency of 8 fim, emission as a function of stellar mass/spectral type, 
we divide the sample into 'early' (BA), 'intermediate' (FG), and late (KM) stars based on 
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their J-K^ colors. The frequency of 8 fim emission in these groups is 6.3% (1/16), < 3.8% 
(0/26), and 0-1.9% (0-3/158). With the exception of the warm debris disk source HD 45435 
(ID 18566), no star plausibly earlier than KO has 8 fim excess. At least one and as many 
as three late-type stars have large 8 fim to 4.5 fim flux ratios indicative of warm dust. If 
two stars (IDs 18601 and 18622) with small log(F[8]/F[4.5]) but red K,-[5.8, 8] colors are 
also classified as excess sources, the frequency of warm dust around late-type stars could be 
as high as 3.2%. Optical spectra is required to confirm the membership status of late-type 
candidate cluster members and thus better constrain the frequency of late-type stars with 
warm dust. 

These results suggest that warm dust surrounding NGC 2232 cluster stars is rare. How- 
ever, the derived frequencies of warm dust for early-type members are based on a small 
sample, which precludes a robust determination of the excess frequency for 25 Myr-old B 
and A stars. Lack of deep x-ray data and incomplete optical spectroscopy prevents robust 
estimates of membership probabilities for intermediate and late-type stars. These data would 
allows us to separate true cluster members from foreground M dwarfs and background red 
giants. 



4-3.2. Frequency of 24 fJ'fn Emission 

To identify sources with 24 fim excess, we consider the sensitivity of the MIPS data 
and the photospheric colors of likely cluster members. For 5a detections with [24] si 10.5, 
MIPS cannot detect the photospheres of cluster stars with spectral types later than ~ G2. 
For early G stars, the Kurucz-Lejeune stellar atmosphere mod els predict K^-[8] ^ -0.020 



and K5-[24] ~ 0.005. Observations of ~ 30 solar-type stars from lRieke et al.l (120081 ) suggest 
K,s-[8] ^ 0.046 and Ks-[24] ~ 0.045. Uncertainties in absolute calibration are small, ~ 1.5% 



(|Rieke et al.l 120081). Older c alibrations suggest K^-L'/M < 0.02 for solar-type stars (e.g. 



Kenyon and Hartmannlll995l ). Thus, we expect small photospheric colors - Ks-[24] < - for 



all cluster stars detected by MIPS. 

Small number statistics and different areal coverage prevent developing an empirical lo- 
cus of photospheric MIPS colors for all spectral types. There are only 38 candidate/confirmed 
cluster members detected by MIPS compared to 200 detected at 8 fim. Therefore, it is dif- 
ficult to define the empirical locus of photospheric MIPS colors from early to late spectral 
types. Because the MIPS coverage is slightly larger than the IRAC coverage at 4.5 fim,, not 
all MlPS-detected stars have measured 4.5fim fluxes. Therefore, using the [4.5] filter as a 
photospheric baseline is not ideal. 



-22- 



Because the expected photospheric colors of B to G stars are small and vary little 
with spectral type, we develop an excess criterion based on the Kg- [24] color and the MIPS 
photometric error: 

{Ks - [24])ots - {K, - [2A])phot > 3(724. (1) 

To account for systematic errors (e.g. uncertainities in the zero-point flux, phase-dependent 
response errors, and uncertainties in the model atmospheres) and the observed dispersion of 
colors, we also require that the excess contribute greater than ~ 15% of the total flux at 24 

{Ks - [24\Us - {Ks - mUot > 0.15. (2) 

From Figure [H it is clear that the frequency of 24 fim excess emission is much higher 
than for 8 fim excess emission. As previously noted, the poorer sensitivity of MIPS results 
in completeness-related bias for cluster members later than ~ G2; F and G stars are also 
typically too faint to have proper motion data and robust x-ray detections. In contrast with 
the A star population, the census of F and G cluster members is likely incomplete. Thus, 
we identify excess sources for FGK stars but do not attempt to derive a frequency of 24 /im 
excess emission for these stars. 

Based on the criteria for excess, 8/15 B and A cluster stars have 24 //m emission from 
circumstellar dust|j. We identify two additional excess sources (IDs 4144 and 5494) with 
J-Kg colors suggestive of FO-KO stars. Other sources with spectral types between FO and 
KO have red K5-[24] colors suggestive of IR excess. Because their colors are less than 3a 
away from the locus of photospheric sources, we do not classify these sources (IDs 11320, 
18591, and 18605) as excess sources. Among the later-type stars, we identify ID 18012 an 
excess source. The spectral types estimated from photometry (~ K7-M0) and spectroscopy 
(MO) are consistent. Thus, ID 18012 is a rare example of a 10-30 Myr-old M star with 
substantial debris emission. Table E] lists properties of the early- type members of NGC 2232 
and late-type members with 24 ^m excess. 



^We count ID 18648 as a late A star from its J-K^ colors, though removing this source from the list of B 
and A stars does not qualitatively change our analysis. 
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4.4. Evidence for a Correlation between Rotation and Debris Disk Emission 
for High-Mass Stars: Massive Primordial Disks May Evolve into Massive 

Debris Disks 



Stellar rotation is another property that can impact the frequency and luminosity of disk 
emission. With studies traditionally focused on near-lR ground-based observations of young 
clusters with primordial disks, the c onnection between disk emission and ste l lar rotation has 



been controversial. Many authors (jStaussun et al.lll999l: iRebuU et al.ll2004j : iMakidon et al. 



20041 ) fail to find a correlation. Others, most notably iHerbst et al.l (120021 ). claim that slower 
rotators are more likely to have disk emission. Because the level of disk emission compared 
to the stellar photosphere is larger in the mid infrared than in the near-IR, Spitzer Space 
Telescope observations of young stars provide a more definitive test of the connection between 
stellar rotation and disk emission. 



Recently, ICieza and Baliberl (120071 ) investigated the stellar rotation-disk emission con- 
nection in two young clusters, NGC 2264 (~ 2-3 Myr) and the Orion Nebula Cluster (~ 1 
Myr). By cor nparing the rotation periods and Sp itzer/IRAC excesses of high-mass cluster 
stars (see also lRebuU et al.ll2006l : llrwin et al.ll2008 ). they identify a clear increase in the disk 
fraction with rotational period. ICieza and Baliben interpret this correlation as evidence in 
favor of the "d i sk locking" model f or rotational evolution of pre-main sequence stars (e.g. 
Shu et al.lll994j : iBouvier et al.lll997l ). where the star loses angular momentum from magnetic 
field interactions with the disk. Thus, the more massive disks (presumably identified by their 
stronger mid-IR emission) are able to spin down the star more substantially. 

If massive primordial disks evolve into massive debris disks, then stellar rotation rates 
should correlate with the level of debris emission. While MIPS data for the (3 Pic Moving 
Group hint at a po ssible connection between the rotation rate and IR excess in debris disks 
( jRebuU et al.ll2008l ). the authors caution that such a correlation is preliminary. The sample 
is dominated by GKM stars that can drive strong winds from their convective atmospheres, 
thus affecting stellar rotation. A more conclusive plan is to compare rotation and MIPS 
excess in B/A star cluster members, which lack convective atmospheres and generally have 
stronger 24 /xm excesses. 

To test this hypothesis, we combine the MIPS data for B/A star cluster members with 
the rotational velocity study of iLevatd (119741 ). The study derives vsini from the He I (4471 
A) and Mg II (4481 A) l ines wit h an uncertainty of ~ 10-20%. We match all but three of 
the members studied by iLevatd . All of these stars are photometrically, spectroscopically, 
and astrometrically confirmed as members. None have been identified as Be stars, though 
the rotation rate of one star, with vsini ~ 400 km s~^, is close to the breakup velocity for 
2-2.5 Mq stars. 
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Figure [12] shows that rotation and 24 fim excess appear correlated. Only one star with 
vsini > 200 km s~^ has a (marginal) excess of ~ 0.27 mag. In contrast, 75% of stars with 
vsini < 200 km s~^ have 24 /im excesses. The MIPS excesses of the slower rotators also 
extend to much larger amplitudes, up to ~ lOx photospheric levels, compared to those of 
the faster rotators (~ 25% of the photospheric levels). The Kolmogorov-Smirnov test yields 
a probability of ~ 1.6% that there is no correlation between rotation and 24 fim excess 
(d=0.8). Both B stars and early A stars have excesses and photospheric emission. The 
mean and median spectral types of stars with vsini < (>) 200 km s~^ are B8.7 and AO 
(B7.6 and B6.7). 

While a larger sample of B and A stars would provide better statistics, these data 
support the picture that massive primordial disks evolve into massive debris disks. However, 
several caveats about rotational velocities for early-type stars complicate our interpretation 
of Figure [121 For example, we only know the component of the rotational velocity along 
our line of sight, so the true rotation rate is not known. While the effects of sini should 
wash out in very large samples, they may be important for the NGC 2232 sample alone. To 
better investigate any connection between rotation and debris emission, these data should 
be combined with data for other similarly-aged A stars to provide a much stronger statistical 
significance to any correlation. 



5. The Evolution of Debris Emission from Planet Formation: NGC 2232 in 

Context 

5.1. Evolution of 8 nm Emission 

The low frequency of 8 fim emission from disks in NGC 2232 is consistent with recent 
Spitzer surveys of other young clusters with substantial populations of warm, terrestrial zone 
debris disks (Table [7]). In h and x Persei (13 Myr), the frequency of 8 fim circumstellar dust 
emission rises from ~ 1.3-2.2% (14/1023 ; 42/1878) fo r B and A stars to > 6% (34/523; 



50/618) for FG stars fICurrie et al.ll2007al ; ICurriell2008[ ). ICurrie et al.l fl2007al ) cite a lower 



limit of ~ 4% for the frequency of 8 nm emission from intermediate-type stars. This lower 
limit is still larger than the upper limit for the early-type excess frequency of B and A stars 
for h and x P^r (r^ 3.2%) and much larger than that for B stars alone (< 1%). Thus, the 
frequency of 8 fim emission is spectral-type dependent at a given age and environment. 

For other young clusters, small number statistics prevent a clear m easure of the va r iation 



of excess 8 fim. emission with spectral type. In NGC 2547 (38 Myr), iGorlova et al.l (120071 ) 



derive excess frequencies of <5% (0/19) for B and A stars, ~ 1.5-5% (1-3/62) for F and G 
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stars, and ~ 1% (3/347) for K and M stars. Several stars studied in iGorlova et al.l have 8 
fim excesses but were not classified as excess sources because they lacked detectable excesses 
at 5.8 fim. This criterion may preclude detection of colder dust where the 8 /im excess is on 
the Wien side of the blackbody. Relaxing the requirement of a clear 5.8 /im excess in NGC 
2547 stars yields a larger frequency of warm dust emission, up to ~ 2% (8/347) for M-type 
stars in this cluster. In NGC 2232, our estimates for the frequency of 8 /im excess - ~ 6% 
(BA), < 4% (FG), and ~ 0-1.9% (KM)- are nearly identical to the results for NGC 2547. 
Within the errors, the low frequency of 8 fim emission among AB and FG stars in NGC 2232 
and NGC 2547 is qualitatively consistent with results for h and x Pg^'- 

Combining results for all four clusters, the frequency of warm dust among 10-30 Myr 
cluster stars at a given age and environment is small. Across all spectral types, the frequency 
of 8 yum emission is < 5-10%. Although NGC 2232 and NGC 2547 have relatively sparse 
populations of stars with G or earlier spectral types (~ 150), the large population in h and 
X Per (> 2000) yields a robust estirnate fo r the frequency of warm dust among early and 
intermediate-type stars (ICurrie et al.ll2007al ). Current data preclude 8 /zm measurements for 
K and M stars in h and x Per- However, the large populations in NGC 2232 (158) and NGC 
2547 (347) demonstrate that warm dust at levels L^/L^ > 10~^ is rare among 10-30 Myr old 
K and M stars. 



5.2. Evolution of 24 nm Emission 



Previous Spitzer studies of 5-100 Myr-old clusters examine the time evolution of 24 
Atm excess emissio n ar ound A stars. U sing somewhat different samples of nearby A stars, 
Rieke et al.l (120051 ) and ISu et al.l (120061 ) show that the maximum amount of ex cess emission 
at a giv en stellar age declines inversely with age from ~ 30-50 Myr to ~ 1 Gyr. ICurrie et al. 



( l2008al ) add results for several young clusters to this sample and demonstrate that the 
maximum level rises for stars with ages of 5-10 Myr, peaks for stars with ages of 10-20 Myr, 
and then declines inversely with time. With new results for NGC 2232 and NGC 2547, we 
can reinvestigate 'the rise and fall of debris disks' and consider, for the first time, the relative 
frequency of debris disks for A stars in young clusters. 

To investigate the evolution of 24 /im emission among A stars in young clusters, we 
consider clusters with good statistics and ages of 5-100 Myr. In addition to our results for 
NGC 2232, we include data from Orion OBla and b (10 and 5 Myr), Upper Centaurus Lupus 
(16 Myr), h and x Persei (13 M yr), NGC 2547 (38 My r ), IC 2602 (50 Myr) IC 2391 (50 
Myr) , and the Pleiades (100 Myr iHernandez et al.l l2006l : ICurrie et al.l l2008al : iGorlova et al. 



2007 



Sieder et al.l l2007l : IGorlova et al.ll2006l ) 
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These data demonstrate a time-dependent frequency of 24 /xm emission from dust among 
young A stars (Figure [T3l) . For stars with ages 5-25 Myr, the frequency appears to rise from 
< 40% for 5-10 Myr old stars to > 50% for 15-25 Myr old stars. For older stars (50-100 
Myr), the 24 /ini excess frequency declines to ~ 10-20%. 

The increasing 24 fim excess frequency from ~ 5 Myr to ~ 10-25 Myr can be explaine d 
within the context of standard models of planet formation (e.g. iKenyon and Bromleyll2008l ). 
During the early stages of planet formation, much of the primordial dust mass may be 
incorporated into > km-sized planetesimals. This growth to larger sizes reduces the disk 
opacity and thus level of emission. Some disks may have undetectable levels of emission. 
Once forming planets reach > 100-1000 km sizes, they stir the leftover planetesimals they 
accrete to high velocities. When planetesimals collide with each other at these velocities, 
they produce copious amounts of dust which is then detectable in the mid infrared. 

Next, we ex ami ne the level of 24 u m excess emission with respect to time (see also 
Rieke et al.ll2005l and lCurrie et al.ll2008al ). focusing on 5-40 Myr-old early-type stars. We in- 
clude the complete sample of B and A stars in NGC 2232. In addition to ou r data, we include 



stars from clusters wi th debris disks: Ori on OBla and b (10 and 5 Myr; iHernandez et al. 



20061), 7] Cha (8 Myr; iGautier et al.ll2008h. the 3 Pic Moving Group (12 Myr; iRebuU et al. 
20oi). h ari d y Persei f^ 13 Mvr: IC urrie et al.l [2008al) . Sco-Cen, and NGC 2547 (38 Myr; 
Rieke et al.ll2005l ). For the Sco-Cen, we include data for Upper Scor pius (5 Myr ) , Uppe r 
Centaurus Lupus (16 Myr), and Lower Centaurus Crux (17 Myr) from lChen et al.l (j2005al ). 
We limit our literature sample to BAF stars. 



Fig ure [Til shows the amplitude of 24 jj,m excess vs. time. In support of ICurrie et al. 
( l2008al ). there is a rise in the level of debris emission from ~ 5 Myr to 10 Myr and a peak 



at ~ 10-20 Myr|. Adding > 30-50 Myr-old clusters JCurrie et al.ll2008al : [Rieke et al.ll2005h 
shows that the emission is consistent with a 1/t decline. The mean level of excess (dashed 
line) for < 20 Myr-old clusters clearly deviates from a t~^ decline (dot-dashed line) observed 
for > 30-50 Myr-old stars. 

These new data improve our understanding of the evolution of debris disk emission at 
~ 20-40 Myr. Data for NGC 2232 suggest a slow decline in debris emission from 10-20 Myr 
to 25 Myr. The mean level of excess emission from the f3 Pic Moving Group (12 Myr) and 
Sco-Cen (16-17 Myr) to NGC 2232 (25 Myr) clearly drops. However, the excess levels for 
NGC 2547 are comparable to those for NGC 2232. The trend in the decline of 24 fim excess 
emission is then slightly shallower than a 1/t (dash-three dots) through ~ 30-40 Myr. As 



^Unpublished da ta for other 5-40 Myr old clusters 
ICurrie et al.ll2008di Currie 2008b, unpublished). 



NGC 2362 and IC 4665 - also support this view 
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shown in iKenyon and Bromleyl (120081 ). debris disk evolution models for A stars predict a 
high level of emission thr ough ^ 30 Myr. Debris enaissioii declines more rapidly for sources 
older than ~ 40-50 Myr jRieke et aljboosi : ISu et alj[2006[ ). 



These results demonstrate that an epoch of detectable, strong mid-IR emission fol- 
lows an epoch with weak or undetectable emission. When observed in an ensemble of stars 
at a range of ages, this sequence manifests itself as an increase of the 24 fim excess fre- 
quency and luminosity with respect to time. This trend is not apparent at Sfim because the 
time for protoplanets to form in the terrestrial region probed by 8/im dust emission is far 
shorter than in the ice giant /Kuiper belt regions which are probed by 24 nm dust emission 
( IKenyon and Bromleyl 120041 ). If 8 fim excess emission is systematically weaker than 24 fim 
excess, scatter in the photometry for intermediate and late-type members may conceal bona 
fide excess sources. 

The high frequencies of 24 /im dust emission from multiple 10-40 Myr-old clusters 
imply that at least the majority of primordial disks pass through the debris disk stage at 
later times. Dust with emission at 24 ^m but not at 8 /xm probes active planet formation 
in regions with equilibrium temperatures of ~ 120 K (jCurrie et al.ll2008al ). which is colder 
than the water-ice condensation temperature in the solar nebula (~ 170 K). 



Within the context of the IKenyon and Bromleyl debris disk evolutionary models, 24 fim 
dust probes active planet formation beyond the ice line, so the frequency of icy planets 
around A stars (~ 1.5-3 Mq at this age) is rji ^ 50-60%. Thus, most ^ 1.5-3 Mq stars 
should form icy planets. 



Discussion 



We have described the first search for debris disks in the 25 Myr-old nearby open cluster, 
NGC 2232. Using ROSAT x-ray observations and optical/IR colors and spectroscopy, we 
identified probable cluster members and cross correlated this list with Spitzer IRAC and 
MIPS data. By comparing these data with results for other 5-40 Myr old clusters, we 
reconstruct the time evolution of 8 ^m and 24 ^m emission from debris disks and probe the 
connection between debris emission and stellar rotation. 

Our analysis identifies one A-type cluster member, HD 45435, with strong warm dust 
emission and 1-3 late- type stars with weakly-emitting warm dust. HD 45435 has dust with 
a temperature of T ~ 350-800 K and a luminosity Ld/ Li, ~ 5 x 10^'^. For the late-type stars, 
the lack of 24 /xm excesses yields good limits for the dust luminosity, L^/ L^, ~ 1 — 10 x 10""^, 
but poor limits on the dust temperature, T < 200-300 K. The vast majority of cluster stars 
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lack evidence for marginal warm dust emission. Combined with results for h and x P^r and 
NGC 2547, our analysis suggests that warm dust emissi on from debris disks a round 10-30 



Myr o l d stars is rare, with a frequency < 10% (see also iMamajek et al.ll2004l : ICurrie et al. 
2007al : ICorlova et al.ll2007l ). 



In contrast, cold dust emission around 10-30 Myr old A- type stars is common. Roughly 
half of all A stars with ages between 10 Myr and 30 Myr have excess emission at 24 /im. 
Upper hmits to the excess emission at 8 /im yield good constraints on the dust temperature, 
T ~ < 200-300 K, and luminosity, Ld/L^ < lO'^. 

These results place interesting limits on the frequency of planet around young stars. 
For most A-type stars with 24 fun excesses, dust emission probes ac tive planet formation 



beyon d the ice line, according to debris disk evolutionary models (e.g. iKenyon and Bromley 
20081 ). Thus, at least 50% of A-type (~ 1.5-3 Mq) stars hkely have icy planets. Warm dust 
emission probes terrestrial planet formation inside the ice line. Thus, < 10% of all stars 
show evidence for rocky planet formation. Despite fairly large samples, it is not yet clear 
whether terrestrial planet formation (a) is rare, (b) produces only weak emission concealed by 
intrinsic photometric scatter, or (c) concl udes much more rapidly tha n icy planet formation, 
as predicted by planet formation models (JKenyon and Bromleyll2004l ). 



For well-sampled clusters with ages 5-50 Myr, the frequency and magnitude of cold dust 
emission among early-type stars depends on stellar age. The frequency of 24 /xm emission 
increases from 5-10 Myr, peaks at 10-30 Myr, and then declines. The typical amount of 
cold dust emission correlates with the frequency of dust emission. Both of the se trends are 
consistent with predictions from debris disk models (JKenyon and Bromleyll2004j ) . Comparing 
the rotation rates of B and A stars with their 24 fim excesses shows that slower rotators 
have stronger excesses. B ecause primordial disk emission is also stronger for slower rotators 
(jCieza and Baliberl 120071 ). there may be an evolutionary link between massive primordial 
disks and massive debris disks. 

Despite our conclusion that icy planets are common around higher-mass stars, we cannot 
infer the masses of planets from the luminosity of the debris disk emission. The objects 
responsible for debris-producing planetesimal collisions have radii > 500 km and masses > 
10~^ M^. While the debris luminosity increases as planets grow to hundreds of kilometers 
in size during runaway growth, the luminosity is independent of the planet size during 
oligarchjo. The observed L^^/L^ depends on the initial disk mass - and thus the total mass 



^Even though planetesimals colhde at higher velocities around more massive planets, the scale height 
of the planetesimals also increases. Thus, the collision rate, which control s the debris luminosity, does not 



increase dramatically (S. Kenyon, unpublished: iKenvon and Bromlevll2008l ) 
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in larger objects - more than the mass of the largest object (JKenyon and Bromleyll2008l ). 



Radial velocity surveys also constrain the frequency of planets around high-mass staro- 
Current surveys show that the frequ ency of planets is h igher around > 1.5 Mq stars than 
around solar and subsolar- mass stars (I Johnson et al.ll2007l ). consistent with theoretical expec- 
tations (JKennedy and Kenyonll2008l ). From the standpoint of debris disk studies, however, 
it is not clear whether a) high-mass stars should make more icy planets than solar-mass 
stars or b) high-mass stars produce more debris than they do around solar-mass stars simply 
because their disks are likely more massive. 

This work highlights two problems that impede our understanding of terrestrial planet 
formation from surveys of young, open clusters with Spitzer observations. The lack of massive 
young clusters prevents robust constraints on terrestrial planet formation around high and 
intermediate-mass stars. The low luminosity of warm debris around low-mass stars limits 
our ability to construct large samples of stars with warm dust. High-precision photometry 
of nearby clusters like NGC 2232 may probe terrestrial planet formation around low-mass 
stars. However, NGC 2232 and most other nearby clusters are not massive enough to provide 
constraints on terrestrial planet formation around intermediate and high-mass stars. The 
only cluster s tudied so far that is clearly massive enough to provide robust constraints is h 
and X Persei ( ICurrie et al.ll2007a| ]bl). 



To investigate the evolution of warm debris emission around high and intermediate-mass 
stars more conclusively, observations of clusters more populous than NGC 2232 are necessary. 
Spitzer Cycle 4 and 5 observations of h and x Persei and other 8-20 Myr old clusters will 
fill this role. Preliminary analysis of MIPS h and x Per data taken during March 2008 show 
a huge increase in the num ber of cluster stars d etected, including several terrestrial zone 
debris disk candidates from ICurrie et al.l (l2007bl ) that previously had MIPS upper limits. 
Observations for h and x Persei will be bracketed by observations of massive 8 and 20 
Myr old clusters, NGC 6871 and NGC 1960, to provide better observational constraints on 
terrestrial planet formation during the critical 10-30 Myr epoch. 

This work justifies the need for follow-up observations of NGC 2232. New optical pho- 
tometry, spectroscopy, and x-ray surveys will provide a far better constraint on cluster mem- 
bership and member properties, which will improve the utility of NGC 223 2 as a laboratory 
for studying debris disk evolution and planet formation. iLyra et al.l (120061 ) showed that op- 
tical colors provide a useful diagnostic of cluster membership; a deep spectroscopic survey 
to identify late-type members is also feasible with multi-object spectrographs like Hectospec 



^"^However, constraining the frequency of icy planets around high-mass stars requires an observing baseUne 
that is much longer than that which is currently available. 
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on the MMT and Hydra at Kitt Peak. Identifying late-type members with a deep Chandra 
survey of NGC 2232 would better constrain the evolution of x-ray activity from 10 Myr to 
50 Myr in addition to identifying late-type, chromospherically active cluster stars. All these 
observations should provide a much more complete census of cluster members, which in turn 
will yield a better analysis of the disk population for late-type stars. 

Finally, the proximity of NGC 2232 to the Sun allows for follow-up ground-based/space- 
based observations of disks, some of which are not possible with distant, populous clusters 
such as h and x Persei, NGC 6871, and NGC 1960. The hot dust emission from ID 18566 
may be resolvable by ground-based interferometers such as the Large Binocular Telescope 
Interferometer. The ESO space telescope Herschel should be able to survey NGC 2232 
for evidence of cold debris disks emitting longwards of ~ 40 /im. With a significant disk 
population established, NGC 2232 should be an excellent target for mid-IR observations 
with the James Webb Space Telescope. 

Although NGC 2232 has been virtually ignored by the star/planet formation commu- 
nity for the past 30 years, this work shows it to be a potentially important laboratory for 
understanding debris disk evolution and planet formation which deserves further study. 

We thank the anonymous referee for a thorough review and helpful suggestions which 
improved the quality of this paper. We also thank Nancy Remage Evans for advice on 
ROSAT archival data and Jesus Hernandez for use of the SPTCLASS spectral-typing code. 
Finally, we thank Ken Rines for taking a spectrum of ID 18012 on short notice. This work is 
supported by Spitzer GO grant 1320379, NASA Astrophysics Theory grant NAG5-13278, and 
NASA grant NNG06GH25G. This paper makes use of the WEBDA open cluster database. 
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Table 1. Photometry Catalogue for Sources with Detected on NGC 2232 field 



a([3.6]) [4.5] <t([4.5]) [5.8] <t([5.8]) [8] 



a ([8]) [24] 



<t([24]) 



1 


96.4706 


-5.3994 


0.00 


13.51 


13.05 


12.86 


0.00 


-99.00 


0.00 


-99.00 


0.00 


-99.00 


0.00 


-99.00 


10.25 


0.15 


2 


96.4708 


-5.4639 


0.00 


11.19 


10.46 


10.19 


0.00 


-99.00 


0.00 


-99.00 


0.00 


-99.00 


0.00 


-99.00 


9.78 


0.14 


3 


96.4790 


-5.1864 


0.00 


15.82 


15.12 


14.90 


0.00 


-99.00 


0.00 


-99.00 


0.00 


-99.00 


0.00 


-99.00 


9.22 


0.06 


4 


96.4816 


-5.3179 


0.00 


10.52 


9.91 


9.71 


0.00 


-99.00 


0.00 


-99.00 


0.00 


-99.00 


0.00 


-99.00 


9.66 


0.09 


5 


96.4836 


-5.0643 


0.00 


9.90 


9.67 


9.56 


0.00 


-99.00 


0.00 


-99.00 


0.00 


-99.00 


0.00 


-99.00 


9.68 


0.09 



Note. — The tabic includes sources with at least one detection in a Spitzer IRAC or MIPS band. Values of 0.00 in the photometry column and 
-99.00 in the photometric uncertainty columns denote sources that were not observed in a given filter. Sources observed in the MIPS band but not 
in the IRAC bands fall outside the IRAC coverage. 



Table 2. Spectroscopy of Selected NGC 2232 sources 



Spectral Typo (T(Spcctral Type) EW{Hq) FWHM(Hq) ROSAT counts ks" 



1050 


96.5720 


-4.5989 


GO. 7 


3.4 


6041 


96.8377 


-4.9424 


G6.3 


2.4 


6639 


96.8691 


-4.8748 


G3.7 


2.0 


10258 


97.0566 


-4.5612 


K6.0 


1.2 


10763 


97.0836 


-4.6977 


G2.5 


2.6 


16505 


97.3680 


-4.6234 


F4.0 


1.4 


18012 


97.4388 


-5.0088 


MO.O 


1.5 


18553 


96.8783 


-4.6588 


F9.8 


1.9 


18554 


96.6436 


-4.5974 


B3.1 


1.3 


18555 


96.8426 


-4.5690 


K1.8 


1.2 


18558 


96.7537 


-4.3556 


B3.9 


1.4 


18559 


96.6810 


-4.7450 


F5.0 


3.0 


18560 


96.5763 


-4.5993 


M3.2 


0.9 


18562 


96.6217 


-4.3049 


M3.6 


1.0 


18566 


96.7940 


-4.7801 


A2.8 


2.5 


18568 


97.0254 


-4.6412 


MO. 8 


0.7 


18570 


96.8788 


-4.5815 


M3.1 


0.8 


18572 


96.8402 


-4.9407 


F7.3 


2.3 


18577 


97.0173 


-4.8984 


K3.9 


1.2 


18578 


96.8713 


-4.8730 


M3.6 


0.9 


18579 


96.9484 


-4.8247 


B7.6 


1.2 


18580 


97.0004 


-4.7930 


M2.9 


1.2 


18581 


96.9899 


-4.7621 


BO. 8 


1.4 


18582 


96.8650 


-4.7564 


MO. 8 


0.8 


18587 


97.1874 


-4.7139 


MO. 5 


0.8 


18590 


96.9365 


-4.7342 


M4.2 


1.1 


18599 


96.9009 


-4.9428 


M3.5 


1.0 


18605 


97.1404 


-4.7829 


G1.4 


2.5 


18606 


96.9688 


-4.7669 


B9.5 


1.3 


18607 


97.0515 


-4.7617 


M3.9 


1.2 


18612 


97.0250 


-4.6057 


M2.4 


0.7 


18613 


97.0584 


-4.5622 


K7.5 


0.8 


18617 


97.0399 


-4.7984 


M2.0 


0.8 


18620 


97.0354 


-4.7400 


M2.2 


0.8 


18621 


96.9177 


-4.7127 


M3.7 


1.2 


18622 


97.0234 


-4.6411 


M3.4 


1.1 


18633 


96.5333 


-4.6282 


B8.8 


1.2 



2.70 


99 


- 


2.40 


99 


- 


0.00 


99 


- 


1.60 


99 


- 


2.00 


99 


- 


4.50 


99 


- 


-0.7 


99 


- 


2.30 


99 


6.1 


5.90 


99 


1.4 


-0.15 


99 


2.1 


5.90 


99 


0.9 


5.00 


99 


1.5 


-7.40 


5.20 


0.9 


-6.90 


5.00 


0.7 


9.90 


99 


1.5 


-1.80 


4.30 


0.9 


-5.10 


5.40 


1.2 


3.90 


99 


2.6 


-0.76 


5.90 


3.6 


-6.30 


6.00 


1.3 


7.20 


99 


2.1 


-3.70 


4.50 


0.8 


4.60 


99 


5.9 


-1.80 


6.00 


1.9 


-2.28 


4.70 


2.2 


-8.80 


5.70 


0.6 


-5.10 


5.30 


0.8 


2.80 


99 


1.0 


11.10 


99 


0.4 


-7.00 


5.00 


0.5 


-5.05 


5.20 


1.2 


-1.50 


6.30 


1.5 


-3.80 


5.10 


1.1 


-9.25 


5.50 


0.6 


-6.06 


4.90 


2.1 


-8.20 


5.80 


2.6 


8.80 


99 


- 



yes 
yes 
no 
yes 
yes 
yes 
yes 
yes 
yes 
yes 
yes 
yes 
yes 
yes 



yes 
yes 
yes 
yes 
yes 
yes 
yes 
yes 
yes 
yes 
yes 
yes 
yes 
yes 
yes 
yes 
yes 
yes 



Note. — FAST spectra of selected stars in the NGC 2232 field. ID 10258 lias a near-IR colors lying just outside the range for candidate members 
but a spectral type possibly consistent with membership. The uncertainties in spectral type (cr(ST)) are determined from the dispersion in spectral 
types computed from each of the spectral indices (see iHernandez et al.ll2004l . for more information). Entries in the FWHM(Ha) column with a 
value of '99' did not have H^ clearly in emission and/or a high enough signal to noise to calculate a reliable full-width half-maximum. 



-36- 



Table 3. Confirmed and Candidate Cluster Members with Spitzer Data 



ID 



RA 



DEC Mem. Type 



150 


96.5238 


-4.7164 


4 


845 


96.5620 


-4.6887 


4 


964 


96.5676 


-4.4612 


4 


985 


96.5687 


-4.4479 


4 


1084 


96.5738 


-4.8108 


4 



Note. — The table includes sources with at least one detection in a Spitzer IRA C or M IPS band. Membership type 1 denotes 
confirmed x-ray active cluster members; type 2 identifies members from Clarja| lll972h . Type 3 lists 'candidate' members as 
determined from V/V-J colors, spectroscopy, and proper motion. Type 4 lists 'candidate' members determined from J/J-Ka 
and J/J-H colors. 



Table 4. Photospheric 2MASS-IRAC/MIPS Colors 



Spectral Typo Ks-[5.) 



Ka-[8] Ka-[24] 



BO 


-0.17 


-0.24 


-0.34 


B3 


-0.11 


-0.16 


-0.21 


BS 


-0.07 


-0.12 


-0.1 


B8 


-0.05 


-0.1 


-0.1 


AO 


-0.02 


-0.05 


-0.1 


A5 


-0.01 


-0.05 


-0.1 


FO 


0.02 


-0.02 


0.005 


F5 


0.02 


-0.02 


0.005 


GO 


0.01 


-0.02 


0.005 


G5 


0.01 


-0.02 


0.005 


KO 


0.01 





0.005 


K5 


0.08 


0.1 


0.1 


MO 


0.21 


0.24 


0.34 


M2 


0.25 


0.30 


0.41 


M5 


0.28 


0.33 


0.41 



Table 5. Sfim Excess Sources in NGC 2232 



Spectral Type log(Frgi /Fr4 51 ) log(Fr8] /Fu 51 )— Photosphere Excess/crrT) 



6540 


- 


K5-M0 


-0.262 


0.163 


9220 


- 


K7-M2 


-0.349 


0.074 


18566 


HD 45435 


A2.S 


-0.168 


0.292 


18613 


- 


K2-K7 


-0.355 


0.074 



9.61 
4.39 
17.24 
4.35 



Note. — Sources with warm, 8 fim excess emission from circumstcllar dust as determined by their 8 f-tm to 4.5 fim flux ratios. Sources with a 
range of spectral types have this range estimated photometrically. Two other sources, IDs 18601 and 18622, have Ks-[8] colors that are greater 
than 3crrgi away from photospheric predictions. However, these sources do not have [8]/[4.5] flux ratios that are more than 3a away from the 
empirically-determined locus of photospheric ratios. 
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Table 6. Early- Type Members and Intermediate/Late- Type Members with 24 fiin Excess 



Spectral Type Ka-[8] Ks-[24] 24 pm Excess? Vsm ! (km s ^) 



18554 


HD 45321 




B3.1 


-0.21 


-0.23 


no 


245 


18558 


HD 45418, 


9 Mon 


B3.9 


-0.23 


-0.21 


no 


230 


18566 


HD 45435 




A2.8 


1.00 


2.24 


yes 


1S5 


18579 


HD 45516 




B7.6 


-0.13 


-0.08 


no 


270 


18581 


HD 45546 




BO. 8 


-0.19 


0.04 


no 


45 


18606 


BD-04 1526B 


B9.5 


-0.05 


- 


no data 


- 


18616 


- 




- 


-0.13 


-0.13 


no 


- 


18630 


HD 45583, 


V682 Mon 


B8 


-0.17 


-0.06 


no 


70 


18631 


HD 45532 




AO 


-0.04 


-0.11 


no 


320 


18633 


HD 45238 




B8.8 


-0.02 


0.27 


yes 


390 


18634 


HD 45399 




AO 


-0.05 


0.80 


yes 


170 


18635 


HD 295102 




B8 


-0.04 


0.53 


yes 


<40 


18636 


HD 45627 




AO 


-0.04 


0.40 


yes 


50 


18638 


HD 45565 




AO 


-0.08 


1.00 


yes 


60 


18639 


HD 45691 




AO 


0.0 


1.50 


yes 


75 


18648 


HD 295100 




A7-F2 


-0.07 


0.37 


yes 


- 


4144 


_ 




G2-G7 


-0.07 


0.39 


yes 


„ 


5494 


- 




F2-P7 


- 


0.75 


yes 


- 


18012 


- 




MO 


0.21 


1.57 


yes 


- 



Note. — Spectral types come from the VIZIER database except those determined in this work. Sources with a range of spectral types have this 
range estimated from their J-K^ colors. ID 4144 is identified as a candidate member from its position on the J/J-H and J/J-Kg color-magnitude 
diagrams. ID 5494 is identified as a candidate member from its position on the V/V-J, J/J-H, and J/J-K^ color-magnitude diagrams. ID 18012 is 
identified as a candidate member from the J/J-H and J/J-K^ diagrams and then confirmed as a member spectroscopically. 



Table 7. Frequency of 8fim Excess from > 10 Myr-old Clusters 



Age Parent Sample Size Freq. BA 



Frcq. FG 



Freq. KM 



h and x Pcrsci 


13 


4,700 


NGC 2232 


25 


238 


NGC 2547 


38 


450 



1.3% (14/1023) 
2.2% (42/1878) 
6.7% (1/15) 
(0/19) 



6.5%(34/523) 
8% (50/618) 
(0/33) 



0-1.9% (0-3/158) 



1.6-4.9% (1-3/62) 1-2.3% (3-8/347) 



Note. — References are D iCurrie et al.l J2007al ) : ICurrid i 20081) ■ 2) this work, and 3) lGQrlQva et al.l J2007l ). Parent sample size refers to the initial 
photometric sample. The sizes of individual subsets are listed in parentheses in columns 4-6. Values for the fraction of excess sources in NGC 
2547 include an estimate on the fraction of interloping intermediate/late-type sources that are not likely to be proper motion members. Values for 
h and x Persei include the sample of sources that are photometrically consistent w^ith cluster membership as well as a subset with spectra. The 
frequency of 8 fim excess around late-type stars for NGC 2232 is highly uncertain due to the incomplete membership list. 



-38- 



4000 



3000 - 



2000 - 



1000 - 




Fig. 1. — (Top) Distribution of IRAC detections with 2MASS counterparts. From right to 
left, the distributions are for the [3.6], [4.5], [5.8], and [8] channels. The source counts peak 
at a magnitude of 16, 16, 15, and 14.5 for the [3.6], [4.5], [5.8], and [8] channels, respectively. 
(Bottom) Distribution of MIPS 24 fim and 70 iim detections with 2MASS counterparts. The 
number counts in the 24 /xm filter peak at [24] ~ 10.5 while those in the 70 /xm filter peak 
at [70] ~ 4. 
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Fig. 2. — Spectra of NGC 2232 cluster members ID 18566 and ID 18662. Both sources have 
very red Ks-[5.8] and K^ - [8] colors. The insert show the Hq, line profile. 
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Fig. 3. — Color-magnitude and color-color diagrams for NGC 2232. Top left: J vs. K^ -[8] 
color-magnitude diagram. Top right: Ks-[5.8]/Ks-[8] color-color diagram. Bottom left: J vs. 
Ks-[24] color-magnitude diagram. Bottom right: K,,-[8]/K.;-[2 4] color-color diagram. T he 
reddening vector is derived from the reddening laws derived by llndebetouw et al.l (l2005l ) 
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Fig. 4. — J/J-Ks color-magnitude diagram identifying probable clust er mernbers i n NGC 
2232. Green triangles denote x-ray quiet cluster members identified bv lClarial (119721): g reen 
squares identify x-ray active sources. Stars rejected as non members by IClarial (119721 ) are 
shown as grey dots surrounded by open diamonds. The sequence of x-ray quiet cluster 
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Fig. 5. — V/V-J color-magnitude diagram showing the positions of x-ray active cluster 
members and bright members compared with other stars. Symbols and lines as in Figure HI 
The cluster members define a clear locus in V/V-J. 
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Fig. 6. — J/ J-H and J/J-K^ color-magnitude diagra n is sho wing positions of confirmed and 
candidate cluster members. Sources from the IClarial (11972) catalog (triangles) and ROSAT 
detections (squares) are 'confirmed' members. Candidate members were selected based on 
V/V-J (asterisks) and J/J-H and J/J-K^ (crosses) color-magnitude diagrams. 
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Fig. 7 . — J-K .. vs. K5-[5.8] and Ks-[8] for NGC 2232 cluster members and field stars in the 



Clarial (Il972l ) catalog. The symbols are the same as in the previous figures. Most members 
shown are 'candidate' members selected from their positions on the J/J-K^ and J/J-H color- 
magnitude diagrams. We also show the J-K^ colors expected for star with a range of spectral 
types (dash-three dots). The maximum photometric errors are shown in the lower right-hand 
corner. Most faint sources have cr[5.8,8] ~ 0.02-0.04. There are 1-2 sources with blue J-K^ 
colors and [5.8] and [8] excess; a second potential population of excess sources is centered on 
J-K, ~ 0.85 and has K,-[5.8, 8] < 0.3-0.4. 
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Fig. 8. — Color-color diagrams for cluster members and IClarial (11972) field stars. The 
symbols are the same as in previous figures, (top left) J-K^ vs. K5-[24] colors: many cluster 
members have 24 fim excesses, (top right) K^-fS] vs. Ks-[24] colors for members and field 
stars. With one exception, the sources with clear MIPS excesses have photospheric IRAC 
colors, (bottom) Ks-[4.5] and Ks-[5.8] vs. Ks-[8] color-color diagrams. Almost all sources 
lie in a locus consistent with photospheric colors from early to late- type stars. The sp ectral 
types shown in the top-left panel (dash-three dots) are de rived from lSiess et al.l (120001 ) using 
the color conversions from iKenyon and HartmannI (119951 ) . 



-46- 



1 








^KMl 


^^^^^H 




„ 


■ _ 


■ 


















^^^^^^^^■^&. 


-, 


V 




r^w 


^^^^^^^^^^^^^^^_Jfl 




^^^E^ 


lii 


i 




a^^H 


r 














^^^^^^^^^^^H^^^K 


















^^^^^^^^^^^^^Bk 




■4* -■" 


B 




1^1 




/^ 






u 










" 


^^^^^^^^^ER 




r 






^H 




■ 








^^^^^^^^^^^^^^^v 














fl 










Ip^^^^ 








^^^1 










^^LlI^ n 








-^ n 














^^^^^^^V* 










^^^^^^^^^U^ 














^^^■^ 






^^^^c^ 












m 


~ " " ■ JH^^llil 


-«■ 






2*^ 










■ " 






Fig. 9. — IRAC images of ID 18566 (left) and ID 18601 (right). Neither sources show 
evidence for source blending that could artifically brighten them in the IRAC bands. 
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Fig. 10.— SED models for ID 18566. (top)Blackbody fits to the SED of ID 18566, using a 
stellar blackbody (dot-dashed line), a hot dust component (dotted line), and a warm dust 
component (dashed line). The blackbody fits reveal a population of hot/warm dust. The disk 
has a fractional luminosity characteristic of luminous debris disks, (bottom) SED modeling 
of the star compared to a transition disk model (dashed-dotted line/crosses), a terrestrial 
zone debris disk model (dashed line/diamonds), and a cold debris disk model (dashed-three 
dots/triangles). The stellar photosphere is shown for reference (dotted line). 
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Fig. 11. — The 8 fim to 4.5 fim flux ratio as a function of J-K^ used to identify NGC 2232 
cluster stars with IRAC excesses. The dash-dotted line shows the linear least-squares fit 
to the flux ratios (rms residuals are ~ 3.92%). Dashed lines identify the 4cr limit beyond 
which we identify excess sources. The dotted lines show the 8a limits. IR excess sources are 
identified as those with flux ratios greater than Aa from the photosphere. Four stars in NGC 
2232 have flux ratios consistent with IRAC excesses from warm dust. From left to right, 
these sources are ID 18566 (J-K, si 0.1), ID 18613 (J-K, th 0.8), ID 6540 (J-K, ^ 0.85), 
and ID 9220 (J-K, ^ 0.9). 
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Fig. 12.— Vszm (km s"^) vs. K,-[24] for B1-A3 stars in NGC 2232. The probability that 
Kg- [24] is uncorrelated with Vsini is ~ 1.6%. Thus, the data show that slower rotators 
typically have stronger IR excesses. 
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Fig. 13. — Frequency of 24 /im excess vs. age for B and A stars in several open clusters 
studied by Spitzer. The frequency of debris emission apparently increases from 5 Myr to ~ 
10-25 Myr before declining after ^ 40 Myr. 
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Fig. 14. — Level of 24 /i?7i excess around BAF stars vs. time for 5-40 Myr old clusters. We 
include Orion Obla and b (10 and 5 Myr; triangles), rj Clia (8 Myr; stars), the j3 Pic Moving 
Group (12 Myr; squares), h and x Persei (13 Myr; diamonds), Sco-Cen (5, 16, and 17 Myr; 
asterisks), N GC 2232 ( 2 5 My r filled circles), and NGC 2547 (38 Myr; small crosses). HR 
4796A, from iLow et al.l (120051 ) is also shown as a cross at t~ 8 Myr and [24]-[24]p/jot ~ 5. 
The data are consistent with a rise in debris emission from 5 Myr to 10 Myr, a peak from 
~ 10 Myr to 20 Myr. For stars older than Ri 40 Myr, t he distribution of 24 /zm luminosity 
defines an envelope which is consistent with a t~^ decline ( iRieke et al.ll2005l ). The dashed line 
shows the mean excess. Predictions for the decline in emission from a steady-state coUisional 
evolution model (dash-three dots) follow a t^^ decline. 



